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ABSTRACT 

We present an analysis of deep multiwavelength data for z « 0.3-3 starburst galaxies selected by 
their 70 /an emission in the Extended- Chandra Deep Field-South and Extended Groth Strip. We 
identify active galactic nuclei (AGNs) in these infrared sources through their X-ray emission and 
quantify the fraction that host an AGN. We find that the fraction depends strongly on both the mid- 
infrared color and rest-frame mid-infrared luminosity of the source, rising to ~ 50-70% at the warmest 
colors (i 7 24/im/i 7 7o ( tim 0.2) and highest mid-infrared luminosities (corresponding to ultraluminous 
infrared galaxies), similar to the trends found locally. Additionally, we find that the AGN fraction 
depends strongly on the star formation rate of the host galaxy (inferred from the observed-frame 70 /im 
luminosity after subtracting the estimated AGN contribution), particularly for more luminous AGNs 
(La. 5_8.o kcV > 10 43 erg s _1 ). At the highest star formation rates (~ 1000 Mq yr -1 ), the fraction 
of galaxies with an X-ray detected AGN rises to « 30%, roughly consistent with that found in high- 
redshift submillimeter galaxies. Assuming that the AGN fraction is driven by the star formation 
rate (rather than stellar mass or redshift, for which our sample is largely degenerate), this result 
implies that the duty cycle of luminous AGN activity increases with the star formation rate of the 
host galaxy: specifically, we find that luminous X-ray detected AGNs are at least ~ 5-10 times more 
common in systems with high star formation rates (> 300 Mq yr ) than in systems with lower star 
formation rates (< 30 Mq yr ). Lastly, we investigate the ratio between the supermassive black hole 
accretion rate (inferred from the AGN X-ray luminosity) and the bulge growth rate of the host galaxy 
(approximated as the star formation rate) and find that, for sources with detected AGNs and star 
formation (and neglecting systems with low star formation rates to which our data are insensitive), 
this ratio in distant starbursts agrees well with that expected from the local scaling relation assuming 
the black holes and bulges grew at the same epoch. These results imply that black holes and bulges 
grow together during periods of vigorous star formation and AGN activity. 
Subject headings: infrared: galaxies — star: formation — galaxies: active - galaxies: starburst 



1. INTRODUCTION 

The observed scaling between the mass of a galaxy's 
bulge and the mass of its central supermassive black 
hole (SMBH) points to a fundamental connection be- 
tween the gr owth of galaxies and their BHs. Recent 
findin g s (e.g.. | Magorrian et al.lll998HFerrarese fc Merritt] 
20001 IGebhardt et all I200fl lAlexander et all I2005bt 
Hopk ins et al.l 120061: Iwild et all I2007D suggest that 
SMBHs and bulges generally grow together; however, 
many of the details are still unclear , particularly at in- 
termediate and high r edshift (e.g., Shields et al.l 120061 : 
I Alexander et all I2008al) . The signatures of SMBH and 

1 Department of Astronomy and Astrophysics, Pennsylvania 
State University, University Park, PA 16802, USA 

2 Sterrewacht Leiden, Leiden University, P.O. Box 9513, 2300 
RA, Leiden, The Netherlands 

3 Department of Physics, Durham University, Durham, DH1 
3LE, UK 

4 Space Science Institute, 4750 Walnut Street, Suite 205, Boul- 
der, Colorado 80301 

5 Pontificia Universidad Catolica de Chile, Departamcnto de 
Astronomi'a y Astroffsica, Casilla 306, Santiago 22, Chile 

6 The Johns Hopkins University, Homewood Campus, Balti- 
more, MD 21218, USA 

7 NASA Goddard Space Flight Centre, Code 662, Greenbelt, 
MD 20771, USA 

8 Department of Physics, Texas A&M University, 4242 TAMU, 
College Station, TX 77843, USA 



bulge growth are both readily observable over a broad 
redshift range, as SMBH growth produces an AGN and 
bulge growth is accompanied by active star formation. 
A simple observable that links these two indicators is 
the AGN fraction as a function of star formation rate. 
A determination of the form of this relation over a 
broad range of redshift would provide new constraints 
on large-sc ale models of galaxy f ormation and evolu- 
tion (e.g.. iDi Matteo et all 120081 : iHopkins et all 120081 : 
lYounger et alJl2009f ). 

The AGN fraction, for a given sample, is the number 
of systems with AGN activity divided by the total num- 
ber of systems in which such activity could have been 
detected (e.g., to some AGN luminosity limit), given the 
sensitivity limits of the observations. The AGN fraction 
provides clues to the duty cycle of SMBH accretion: a 
higher fraction implies that the SMBHs spend less time 
in inactive states relative to that spent in active accret- 
ing states. Therefore, any dependence of the AGN frac- 
tion on star formation rate would imply that this duty 
cycle is related to the intensity of star formation. In par- 
ticular, studies that identify AGNs using optical spec- 
tra have shown that the luminous AGN fraction in all 
galaxies at z ~ is on the order of 5-15% percent (e.g., 
iKauffmann et all 120031 : iFrancis et al.l I2004D , whereas in 
massively star forming galaxies, such as the sub-mm 
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galaxies (SMGs) st udied by I Alexander et all (|2005bO and 
lLaird et all (|20100 at z « 2, the fraction is estimated 
to be considerably higher: Alexander et al. find a frac- 
tion of 38±io% using X-ray and radio data, and Laird et 
al. derive a somewhat lower fraction of (20-29) ±7% us- 
ing X-ray-selected AGNs. Although these numbers agree 
within errors, their factor of ~ 2 difference points to the 
need for additional studies of the AGN fraction in this 
high SFR regime, using larger samples and different ap- 
proaches. Additionally, the detailed form of the AGN 
fraction between the two extremes of star formation rate 
is currently poorly constrained. 

Due to dust that absorbs the UV emission from young 
stars and re-emits it at long wavelengths, a galaxy's mid- 
to-far-infrared emission is commonly used as a tracer 
of its star formation activity^ Recently, very deep 
Spitzer MIPS data have become available for the deep- 
est Chandra X-ray fields, which together provide sen- 
sitive X-ray and mid-to-far-infrared observations that 
are ideal for identifying AGN activity to luminosities of 
Lbo\ ~ 10 42 erg s _1 and dust-obscured star formation to 
star-formation rates of ~ 10 M© yr _1 at z ~ 0.5. These 
data allow one to trace luminous star formation and AGN 
activity in the distant universe and to investigate how the 
AGN fraction depends on the star formation rate. 

However, it is well known that AGNs are associated 
with dusty "tori" that are often luminous in the mid- 
infrared. Therefore, AGN may also contribute signifi- 
cantly to the total infrared emission when present. A 
number of studies have investigated the contribution 
from AGN-powered emission to the infrared flux in lu- 
minous infrared sources. Such sources are generally di- 
vided into two subclasses by their integrated 8-1000 /im 
luminosity (denoted £ir): luminous infrared galaxies 
(LIRGs, 10 11 < Lir < 10 12 Lq) and ultraluminous 
infrared galaxies (ULIRGs, Lir > 10 12 Lq). Among 
the general population of luminous infrared sources, star 
formation appears to be the dominant power source of 
the mid-to-far-infrared emission in most objects. For 
example, using diagnostics based on mid-infrared emis- 
sion lines and the stre ngth of the 7. 7 £tm PAH feature 
of z < 0.15 ULIRGS, iGenzel et all (fl998h found that 
star formation likely powers most of the 8-1000 /im lu- 
minosity of 12 of the 15 ULIRGs they studied . Us - 
ing similar mid-infrared diagnostics, iHouck et al.l (|2007f ) 
found that for the majority of their sample, selected by 
24 /mi flux and mostly at z < 0.25, the mid-infrared 
luminosity is dominated by star formation. In a sam- 
ple of 43 0.1 < z < 1.2 objects selected by 70 /zm flux, 
ISvmeonidis et all (2008) fit a variety of starburst and 
AGN-powered emission models to IRAC and MIPS pho- 
tometry and IRS spectra and found that all but one ob- 
ject in their sample are starburst dominated. 

In a study of high-luminosity systems at z < 0.5, 
iTran et all (|2001l) found that at luminosities below Lir ~ 
lO 12 5 Lq, starbursts (identified in ISO spectra by their 
strong mid-infrared PAH emission) are the dominant 
power source in local ULIRGs, but at higher luminosi- 
ties, AGN are often the dominant emission source. How- 
ever, in a study of Spitzer IRS spectra of a sample of 107 

9 In this paper, we define mid-infrared emission to be emission 
at rest-frame wavelengths < 40 /im, and far-infrared emission to 
be at rest-frame wavelengths > 40 fim. 



ULIRGs, iDesai et al.l (|200l found that even the most 
luminous high-redshift ULIRGs often have strong PAH 
emission, indicative of large starbursts th at may be ab- 
sent l ocally. Lastly, the recent study of iVeilleux et al.1 
(2009), which also used IRS spectra of ULIRGs to esti- 
mate the relative contributions of AGNs and starbursts, 
found that the average AGN contribution to the bolo- 
metric (not far-infrared) luminosity of local {z ~ 0.3) 
ULIRGs is ~ 35-40%. However, among far-infrared 
sources with the most luminous AGNs, namely quasars, 
AGN-powered e mission can dominate. For example, 
iShi et all (|2007| ) used the mid-infrared PAH emission 
to infer SFRs in three samples of AGNs: PG quasars, 
2MASS quasars, and 3CR radio-loud AGNs. They found 
that the average contribution of star formation to the 
70 /im emission ranges from 25%-50%, depending on the 
AGN sample. Therefore, one must be careful to account 
for the AGN contribution to the infrared when using it 
to derive star formation rates. 

In this paper, we investigate the growth of SMBHs 
and their host galaxies in a complete sample of starburst 
galaxies, constructed from fields with extremely deep 
multiwavclcngth coverage, and determine the X-ray- 
detected AGN fraction across a broad range of star for- 
mation rate and redshift. Briefly, our sample is con- 
structed using the following approach: 

1. Since mid-to-far- infrared observations sample the 
bulk of reprocessed emission from young stars, we 
use deep mid-infrared (70 /im) data to construct a 
representative sample of star-forming galaxies. 

2. We use deep X-ray observations to identify AGNs 
above a given X-ray luminosity. 

3. In such sources, we use a variety of empirical AGN 
SEDs, scaled by the AGN bolometric luminosity 
estimated from the X-ray emission, to estimate the 
AGN contribution to the infrared luminosity. 

4. Lastly, we use the net infrared luminosity, corrected 
for AGN-powered emission, to estimate star forma- 
tion rates. 

Using this sample, we calculate the X-ray-detected AGN 
fraction above a given limiting X-ray luminosity as a 
function of the SFR, mid-infrared color (a proxy for 
dust temperature), and mid-infrared luminosity, and we 
examine the relative growth rates of the galaxies and 
their SMBHs in distant starbursts. The following sec- 
tions describe in detail each of these steps. We adopt 
Ho = 70 km s- 1 Mpc" 1 , fi A = 0.7, and fl M = 0.3 
throughout. 

2. SAMPLE CONSTRUCTION AND PROPERTIES 

2.1. Mid-Infrared Data 

Samples were drawn from two fields with deep X-ray 
through infrared coverage: the Extended- Chandra Deep 
Field-South (E-CDF-S), which includes the ~ 2 Ms 
Chandra Deep Field-South (CDF-S), and the Extended 
Groth Strip (EGS). The primary sample of star- forming 
galaxies and AGNs was constructed using all sources with 
Spitzer MIPS detections at 70 /im in the Far-Infrared 
Deep Extragalactic Legacy (FIDEL) survey. The FI- 
DEL data comprise very deep coverage of > 90% of the 
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E-CDF-S and EGS fields at 24 fim and 70 fim. Source 
catalogs were created from the DR2 mosaic imaged us- 
ing the DAOPHOT toolQ Aperture fluxes measured by 
DAOPHOT wer e corrected u s ing th e point source func- 
tion derived by iFraver et al.l (|2006f ) to derive the total 
fluxes. No color corrections were performed, as they are 
expected to be < 10% for the bulk of our sources. We 
use the FIDEL 70 fim catalog as the basis of our sam- 
ple because emission at this wavelength should suffer less 
from spectral complexity and have a smaller AGN con- 
tribution for sources at redshifts up to ~ 3 than emission 
at 24 /im. For example, at z > 1, 24 /jm observations 
would sample rest- frame emission at A < 12 /im where 
complex spectral features from PAH emission and silicate 
absorption ar e present in the spectra of many LIRG S and 
AGNs (e.g.. IWeedman et all 120051: lArmus et all 120071: 
IDesai et al.ll2007D . 

Before attempting to identify the counterparts at other 
wavelengths, we made a cut at a 70 /im signal-to-noise ra- 
tio (S/N) of 3, above which the median positional uncer- 
tainty in the 70 /im sources is < 3" (although it reaches 
~ 8" near our adopted S/N cutoff) and completeness is 
high (ps 85%). The S/N and completeness of the 70 fim 
detections were estimated using simulations in which fake 
sources were inserted in the images and their fluxes re- 
covered. To exclude regions of very low exposure near 
the survey edges, where spurious sources are more com- 
mon even at S/N > 3, we also made a cut at an exposure 
time of 1000 s (exposure times were taken from the ex- 
posure maps provided with the FIDEL catalogs). With 
this cut on exposure time, the total areal coverage of re- 
gions with both deep X-ray and mid-to-far-infrared data 
is ps 1100 arcmin 2 in the E-CDF-S and ps 1400 arcmin 2 
in the EGS. The S/N and exposure-time cuts result in 
total source numbers of 567 and 725 in the E-CDF-S and 
EGS, respectively. The 70 /im flux limit for our sample 
varies with the exposure time (by a factor of up to ~ 5) 
and reaches a minimum flux density of ps 1.8 mJy in the 
GOODS-S region of the E-CDF-S. 

2.2. Source Cross Matching 

Due to the large point spread function of the Spitzer 
MIPS instrument and to source blending, 70 /im source 
positions can be uncertain by large amounts (the sim- 
ulations described above give errors of up to ~ 8" 
for a source with S/N = 3). To minimize the num- 
ber of spurious counterparts, we performed a cross 
match between the 70 /im and 24 /im catalogs using 
a probabilistic m atching method (described in detail in 
ILuo et all 120101 ) that takes into account both the es- 
timated source positional errors and their fluxe s (e.g., 
ISutherland k Saunders! fl99l ICilieei etail l200l . This 
method tends to recover a larger fraction of true counter- 
parts than the standard method that uses a single fixed 
matching radius and ignores flux information when se- 
lecting among possible counterparts. Using this method, 
we find that the expected fraction of spurious 24 fim- 
to-70 /im cross matches is < 10%, a value ps 40% lower 
than that obtained with the standard method using a 
fixed matching radius that recovers approximately the 

10 See |http : //irsa. ipac . ca ltech . edu/data/SPITZER/FIDEL 
See http://www.star.bris.ac.uk/~mbt/daophot 



same number of total counterparts P^l With this method, 
527 E-CDF-S and 678 EGS 70 fim sources with S/N > 3 
matched to a 24 fim source. Since the 24 fim data are 
much deeper than the 70 fim data for these fields (by a 
factor of > 30), we can reasonably expect that all 70 fim 
sources should have 24 fim counterparts. Indeed, a vi- 
sual inspection of the ps 6% of 70 fim sources that lack 
24 fim matches suggests that most suffer from signifi- 
cant blending that has distorted their shapes (and hence 
centroids). By requiring a 24 fim match for each 70 fim 
source, we eliminate these problematic sources that likely 
have positions and fluxes in error by large amounts. 

Additionally, ps 25% of the 70 fim sources with identi- 
fied 24 fim counterparts have more than one 24 fim source 
within ps 4". Such multiple matches could imply signif- 
icant blending is present in the 70 fim images, possibly 
leading to spurious cross matches and to misestimates of 
the 70 fim flux. To minimize spurious cross matches, we 
can use the observed distribution of 24/70 fim color for 
sources with unambiguous counterparts to select reliably 
the correct counterpart from multiple matches. To this 
end, we compared the 24/70 fim colors of each possi- 
ble counterpart to the observed distribution of colors for 
70 fim sources with single matches. If two counterparts 
had colors that differed significantly (by > 2a) from the 
mean and the reliabilities determined from the proba- 
bilistic matching process were similar for the two, the 
counterpart with the color closer to the mean was cho- 
sen. In this way, we chose a different counterpart from 
that preferred by the probabilistic matching in ps 5-10% 
of cases (depending on the field) with multiple matches. 
Additionally, in some cases with multiple nearby 24 fim 
sources, blending in the 70 fim images can be significant. 
However, we find that on average only ps 5% of the 24 fim 
sources have 70 fim counterparts. Given that we find 
multiple 24 fim sources for ps 25% of the 70 fim sources, 
each with an average of 2.2 24 fim sources within ps 4", 
we expect blending to have a significant effect on the 
measured 70 fim flux in only 1.2 x 0.25 x 0.05 ps 1.5% of 
the 70 fim sources. At this level, the presence of blended 
sources should not affect our results significantly, and we 
therefore do not attempt to correct for them further. 

Using this 70 fim sample as the basis, we cross matched 
the sources with the following X-ray, optical, and in- 
frared catalogs. For th e E-CDF-S , we u sed the Chan- 
dra X-ray ca t alogs of ILuo et al.l (|2008l CDF-S) and 
iLehmer et all (|2005l E-CDF-S); the COMBO-17 (Clas- 
sifyin g Objects by Me dium-Band Observations in 17 fil- 
ters; IWolf et alJ l 200"l. MUSYC (M ultiwavelength Sur- 
yey by Yale-Chile; IGawiser et all f2006h. and MUSIC 
(|Grazian et al.l [20061) optical and near-infrared catalogs; 
and the SIMPLE {Spitzer IRAC/MUSYC Public Legacy 
survey in the Ext e nded- CDF-S) mid-infrared catalog 
(jDamen et al J 12009( 1. For the EGS we us ed the Chan- 
dra X-ray catalogs of lLaird et ail (|2009H . the optical 

12 The number of expected spurious matches was estimated by 
shifting one catalog relative to the other by 15—60" in 100 differ- 
ent directions, each time cross correlating the catalog sources. The 
average of the resulting number of matches was taken to be the 
expected number of spurious sources for a given search radius. We 
note however that this method likely overestimates the fraction of 
false matches, as we find that the typical distance to a counter- 
part is much less than the search radius when the cross match is 
performed at the real positions. 
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TABLE 1 

Number of sources and counterparts. 



Field 



70 /im 24 /im Optical X-ray Near-IR IRAC 



E-CDF-S 
EGS 



564 
725 



527 
678 



449 
573 



81 
77 



386 
483 



442 
475 



CFHTL33 and CFH12K (Coil ct al. 2004) catalogs, the 
near-infrared Palom ar catalog (iBundv etaL|[2006[ ). and 
the IRAC catalog of iBarmbv et al.l (|2008[ ) 

Since all of our 70 /im sources are required to have 
counterparts in the 24 /im data, we can use the more pre- 
cise 24 /im positions (with typical uncertainties of ~ 1") 
when matching to the optical catalogs, thereby reducing 
the number of likely spurious optical-to-70 /jm matches. 
Additionally, since many of our sources should have ra- 
dio counterparts (since radio emission is often associated 
with starbursts and AGNs), we can further refine the po- 
sitions using the highly accurate astrometry of the radio- 
source catalogs available from VLA surveys of the fields 
(with typical source positional uncertainties of < 0'.'2). 
To this end, we performed cross matching (again using 
the probabilistic matching metho d) between the 24 /im 
sources and the 20 cm catalogs of iMiller et al.l (|2008| 1 in 
the E-CDF-S and version 1.0 of the AEGIS20 catalog 
in the EGS. We found radio matches to 219 (» 41%) 
E-CDF-S sources and 166 (« 25%) EGS sources (note 
that the 20 cm EGS survey does not cover the entire 
field). Due to the low spatial density of radio sources in 
these surveys, we expect < 1% of these matches to be 
spurious. 

We next searched for optical counterparts to the 70 /im 
sources by 24 /im or radio source position. For 70 /im 
sources with an identified radio counterpart, we used a 
matching radius of 0'.'5 (due to the much smaller po- 
sitional uncertainties, probabilistic matching was not 
used). For sources without a radio counterpart, we used 
a radius of 1". When multiple optical matches occur, we 
selected the source with the smallest separation. With 
these matching radii, we expect < 7% of the optical-to- 
70 /im matches to be spurious. We then cross matched 
the X-ray source positions to the 24 /im or radio source 
position, using a variable matching radius that depends 
on the positional uncertain ty of the X-ray s ource. We 
followed the method used in lLuo et al.l (|2008f ) and define 
the matching radius as r = 1.5Ax, where Ax is the X-ray 
positional uncertainty given in the catalogs. Lastly, we 
used a matching radius of 0'.'75 to identify matches to 
mid-infrared Spitzer IRAC sources and a radius of 1" to 
match to UV GALEX sources. Table Q] summarizes the 
results of the cross matching. In total, the final sample 
comprises 1022 unique 70 /im sources with identified op- 
tical counterparts, of which 158 have an identified X-ray 
counterpart. 

2.3. Redshift Estimates 

To obtain redshift estimates for our 70 /im sources, 
we followed the process described above to cross 
match (using a 0'/5 radius) the optical source positions 
to all publicly available spectroscopic redshift cata- 
logs of the E-CDF-S and EGS (e.g., ICristiani et all 

13 See http://www.cfht.hawaii.edu/Science/CFHLS/ 

14 See http://www.roe.ac.uk/~rji/aegis20/ 



2000; iBunker et all 120031: iDickinson et al] [200 



Le Fevre et al.ll2004t iStanway et al j|2004t i Strolger et al] 
20041: Ivan der Wei et al.ll2004t iDavis et alj|2007l) . When 
more than one spectroscopic redshift was available for a 
given source, the redshift of the higher quality (judged 
by the quality flags provided with the catalog) was used. 
If two spectroscopic redshifts were deemed of equal 
quality or the quality was unknown, the average was 
taken. In such cases, the difference between the two 
redshifts was typically < 10%. In total for all three 
fields, 408 (» 40%) of 1022 sources have high-quality 
spectroscopic redshifts determined from two or more 
spectral features. 

Although the majority (~ 60%) of our sources lack 
high-quality spectroscopic redshifts, almost all have high- 
quality photometric data in multiple bands from near- 
UV to mid-infrared wavelengths. Although Wol f et al.l 
(2004) produced a high-quality photometric redshift cat- 
alog for the entire E-CDF-S, new near- and mid-infrared 
data (from the MUSYC JHK and SIMPLE IRAC sur- 
veys; see Table[T]for the fraction of 70 /im sources in each 
field with near- and mid-infrared detections) have re- 
cently become available that should be particularly help- 
ful in deriving accurate photometric redshifts for sources 
at higher redshifts (z > 1.4), which are generally not 
available in the Wolf et al. catalog^ Additionally, no 
photometric redshift catalog is currently publicly avail- 
able for the EGS. Therefore, we can use all existing pho- 
tometric data to obtain additional and improved photo- 
metric redshifts for our sample. To this end, we used the 
publicly available Z urich Extragalactic Bay esian Redshift 
Analyzer (ZEBRA; iFeldmann et al.l l2006) to derive red- 
shifts for the sources that lack a spectroscopic redshift, 
except for CDF-S X-ray sources, for whic h we use the 
photometric redshifts of ILuo et ail (|2010() , derived us- 
ing a very similar method but with a more sophisticated 
treatment of the photometry (e.g., including upper lim- 
its for non-detections). Details of the parameters used 
as input to ZEBRA to derive the photometric redshifts 
and estimates of the quality of the resulting redshifts are 
given in the Appendix. 

In general, the most secure photometric redshifts are 
those for bright sources (tur < 24 AB mags), which com- 
prise the great majority (« 93%) of our final sample of 
1022 sources. In Figure [I] we compare the photometric 
redshifts derived by ZEBRA to spectroscopic ones for the 
subsample of 408 sources in our final sample with spec- 
troscopic redshifts. Assuming the spectroscopic redshifts 
are accurate and the spectroscopic sources are represen- 
tative of the entire sample, ~ 97% of our sources will 
have redshifts with \ztme — z\/(l + Ztruo) < 0.2. For 
AGNs only (identified following fegj) . of which 55 of 108 
have spectroscopic redshifts, the fraction of such sources 
is ss 94%. However, the spectroscopic sample is unlikely 
to be fully representative, and we therefore expect that 
the true fraction of sources with incorrect redshifts will 
be higher by roughly a factor of three (see the Appendix 
for details), but still within acceptable limits (< 10%). 
Lastly, it should be noted that an important conse- 



15 Recently, Cardamone ct al. (2010) used these and other data 
to derive high quality photometric redshifts for the E-CDF-S. We 
have verified that our results do not change when their redshifts 
are used instead. 
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Fig. 1. — Histograms of Az/(1 + z spC c), where Az = z spC c — z p hot 
for all 70 lira sources with spectroscopic redshifts (open region) 
and for AGNs only (filled region). Vertical lines denote Az/(1 + 
z spcc ) = ±0.1 and Az/(1 + z spcc ) = ±0.2. 

quence of obtaining redshift estimates for all of the 70 /im 
sources in our sample with identified optical counterparts 
is that our study will not suffer from strong biases related 
to incomplete redshift coverage. For example, an incom- 
plete sample based on redshift survey s that target AGNs 
or X-ray sources preferentially (e.g.. IZheng et all I2004T) 
would result in an artificially high AGN fraction. 

2.4. Mid- Infrared Luminosities and Colors 

The rest-frame mid-infrared luminosity of a typical 
starburst galaxy or AGN is dominated by reprocessed 
emission from dust. Such emission gives a direct mea- 
sure of the strength of the star formation or AGN emis- 
sion that the dust reprocesses. Therefore, it is of interest 
to investigate how the AGN fraction relates to this lu- 
minosity. To derive the rest-frame luminosities, we first 
constructed observed-frame mid-infrared SEDs from the 
available Spitzer data, which span observed-frame wave- 
lengths from 3.6 /jm to 70 /zm (due to the large posi- 
tional uncertainties inherent to the 160 data, these 
data were not used). The observed SED was then shifted 
to the rest frame of the source using its redshift. We 
then used linear interpolation in log space to derive the 
monochromatic luminosity at a rest-frame wavelength of 
30 /j,m (L30 ; model SED were not used, as a variety of 
AGN, starburst, and hybrid sources that are difficult to 
model are expected to be present). This wavelength was 
chosen to lie within the wavelength coverage of the ob- 
served SEDs of most objects, negating the need for large 
extrapolations. We show the rest-frame mid-infrared lu- 
minosities as a function of redshift in Figure [5^. It is 
clear from this figure that the sensitivity limits of the 
FIDEL survey are such that our sample is roughly com- 
plete only for sources with L30 > 10 12 L@. Below this 
luminosity, the completeness varies with redshift, being 
w 100% at L 3 o > 10 11 L Q to z = 1 and at L 30 > 10 10 L Q 
to z = 0.5. 

Additionally, the mid-infrared color (typically calcu- 
lated using the observed-frame fluxes in the IRAS bands 
as f 25^111/^60^111) of a source gives an indication of the 
temperature of the emitting dust: higher ratios indicate 
relatively more emission at shorter wavelengths, indica- 
tive of emission from warmer dust. Cooler dust tem- 
peratures are likely to be associated with dust heated 
by young stars, whereas warmer temperatures are more 
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Fig. 2. — The distributions of rest-frame 30 fira luminosity (a), 
mid-infrared color (b), and star formation rate (c) of th e 70 /im 
sample as a function of redshift. AGNs (selected in i|2.6|l are indi- 
cated by diamonds. AGNs with net 70 lira flux densities that fall 
below S/N = 3 (after the contribution from the AGN is subtracted) 
have been removed from the lower panel (a total of 21 sources; see 
[|33] for details). 

likely to be indicative of dust heated by AGN emission 
(e.g.. Ide Griip et alJfl985l ISanders et at|[l988l) . There- 
fore, it is of interest to examine the mid-infrared col- 
ors for our sample. In Figure [2]b, we plot the ratio of 
24/70 /im flux (which, for our purposes, we consider to 
be equivalent to the ratio of 25/60 /im flux) against the 
redshift, and in Figure [3J we plot it against the 70 /jm 
flux. For local sources, ratios below log (^24^111/^70^™) ~ 
—0.7 are generally indicative of emission from cool dust, 
whereas higher ratios are indicative of w arm dust (e.g., 
Ide Griip et ; aHll985l ISanders et al.lfl988h . 

It should be noted that the observed mid-infrared color 
for a given galaxy is a strong function of the redshift: at 
higher redshifts (z > 1.5), the portion of the spectrum 
measured by observed-frame 24 /im emission suffers from 
increasing spectral complexity due to the possible pres- 
ence of strong absorption and emission features below 
rest-frame wavelengths of rs 10 /im. Additionally, there 
is evidence that high-redshift (z > 1.5) galaxies exhibit 
stronger PAH emissio n features than local galaxies of the 
same luminosity (e.g.. IMurphv et al]|2009f) that will add 
further redshift-dependent changes to the color. There- 
fore, it is difficult to infer directly a dust temperature 
from the mid-infrared color in sources at z > 1.5 (e.g., 
iPapovich et "all [20071: IMullanev et al.1 [20091) . In gOJ we 
investigate how the inclusion or exclusion of AGN hosts 
with warm mid-infrared colors affects our results. 
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Fig. 3. — The ratio of 24/70 /im flux versus the 70 /im flux. 
AGNs are indicated by diamonds. 

2.5. X-ray Properties 

The purpose of this study is to quantify the AGN frac- 
tion in mid-infrared sources; therefore, a reliable means 
of identifying the bulk of the AGN population is crit- 
ical. Since AGNs are one of only two types of lumi- 
nous X-ray point sources in the distant universe (the 
other being starburst galaxies), and X-rays are not read- 
ily absorbed by surrounding material, the X-ray ob- 
serva tions are extremely effici ent at identifying AGNs 
(e.g.. iBrandt fc Hasingerll2005[ ). In particular, the rest- 
frame 0.5-8.0 keV luminosity and the hard-to-soft X-ray 
band ratio are useful properties in distinguishing between 
AGNs and starbursts. We therefore use the Chandra 
X-ray source catalogs of the CDF-S, E-CDF-S, and EGS 
to calculate the X-ray properties of our sources. 

For the E-CDF-S fi eld, we u s e the 2 Ms CDF-S 
X-ray source catalog of ILuo et all (12008ft and the 250 k s 
E-CDF-S X-ray source catalog of iLehmer etall (|2005l ). 
For t he » 200 ks EGS. we use the X-ray source cata- 
log of ILaird et all (|2009f) . Since a number of differences 
exist between the EGS catalog and the other two cata- 
logs in the band definitions used to measure counts and 
fluxes (e.g., the full band is defined as 0.5-7.0 keV in the 
EGS catalog and 0.5-8.0 keV in the CDF-S and E-CDF-S 
catalogs), we used a simple power-law model to convert 
counts and fluxes to a uniform system. For these conver- 
sions, we first derive the effective power-law index from 
the band ratio given in the EGS cata log following the 
method used in §3.3 of lLuo et all (|2008[ ). Briefly, we find 
the power-law model (including an assumed Galactic col- 
umn density) that reproduces the observed band ratio. 
For sources with a l ow number of c ounts (< 30 counts 
total; for details, see ILuo et al.ll2008T ). we adopt T = 1.4, 
a value representative of faint sources. We then use the 
effective power-law index to convert counts measured in 
the EGS bandpass to the corresponding E-CDF-S band- 
pass. X-ray luminosities and band ratios (the ratio of flux 
in the 2.0-8.0 keV band to that in the 0.5-2.0 keV band) 
were then calculated directly from the catalog fluxes. 

2.6. AGN Identification 

To identify AGNs among the X- ray sources, we fo l- 
low the identification criteria used by Bauer "eTal1(l200l . 
which we outline briefly here. AGNs were first identified 
based on their intrinsic, rest-frame 0.5-8.0 keV luminosi- 
ties. An estimate of the intrinsic absorption is needed 



to derive this luminosity. By assuming that the AGN 
X-ray spectra are well represented by an intrinsic power- 
law with a photon index of 1.8, we can use the band ratio 
(the ratio of counts in the 2-8 keV band to the 0.5-2 keV 
band) to derive a basic estimate of the intrinsic Ah (see 
i i3.II for details of the fitting procedure) . Sources with 
rest-frame Lq. 5-8.0 keV ^3x I0 42 erg s _1 are likely to 
be AGNs, since starbursts generally have luminosities be- 
low Lo. 5-8.0 keV ^ 10 42 erg s _1 . However, to ensure that 
luminous starbursts are not misclassified as AGNs, we 
calculated the predicted 2-10 keV luminosities from sta r 
formation from the relations of lPersic fc Rephaelil (|2007l) . 
Persic et al. find that the scaling relation for ULIRGs is 
different than that of lower-S FR objects. Therefore, we 
use the following relation from lPersic fc Rephaelil (|2007l ) 
to estimate the 2-10 keV luminosity due to star forma- 
tion for systems with SFR < 100 Mq yr -1 : 

L(2 - 10 keV) = 3.8 x I0 39 SFR _ 1 erg s" 1 (f ) 

M Q yr 

For systems w ith SFR > 100 M Q yr _1 , 

iPersic &: RephaeTl ([2007ft found a somewhat differ- 
ent scaling was preferred: 

L(2 - 10 keV) = 0.75 x 10 39 SFR , erg s" 1 (2) 

M e yr 

For the purposes of this calculation, we determined the 
star formation rate following §3.4l bv assuming, conserva- 
tively, that the entire observed 70 (im flux is due to star 
formation. The predicted rest-frame 2-10 keV luminos- 
ity was converted to an observed-frame 0.5-8 keV flux 
assuming a T — 2 power-law spectrum and the source 
redshift. We then classified as AGNs all sources with 
both Lq, 5_8.o keV ^ 3 x 10 42 erg s _1 and with an ob- 
served luminosity > 3 times that predicted from the star 
formation rate. 

Next, sources were classified by the hard-to-soft X-ray 
band ratio: sources with band ratios above 0.8 (corre- 
sponding to effective photon indices T < I) were clas- 
sified as AGN, as starbursts almost always have softer 
X-ray spectra (with T > 1). Lastly, in addition to these 
purely X-ray-based criteria, we also use the X-ray-to- 
optical flux ratio as a further discriminator of AGN ac- 
tivity. We classify sources with /o. 5-8.0 kcv/fa > 0.1 as 
AGNs, where fa is the i?-band flux. Using these criteria, 
we identified AGNs in 108 (« 10%) of the 1022 70 /im 
sources for the combined E-CDF-S and EGS fields. The 
majority of identified AGNs meet more than one selec- 
tion criterion. In particular, of 108 AGNs, 12 were identi- 
fied uniquely using Lo. 5-8.0 keV 3 x 10 42 erg s _1 (which 
identified 94 AGNs in total), 9 using T < 1 (which identi- 
fied 47 AGNs in total), and 7 using /o. 5-8.0 kcv/fn > 0.1 
(which identified 83 AGNs in total). Therefore, we ex- 
pect that few, if any, non-AGNs arc misidcntified as 
AGNs in our sample. 

The remaining 50 non-AGN X-ray sources have prop- 
erties consistent with starbursts and are considered to 
be such in the following analysis. The fraction of X-ray 
sources identified as starbursts in our sample (« 30%) 
is higher than typical starburst f ractions determine d for 
flux-limited X-ray samples (e.g.. iBauer et~aT1l2004l . find 
w 10-20% of the CDF-S X-ray sources are starbursts). 
However, these samples were selected differently from our 
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mid-infrared-selected sample, which naturally includes a 
high fraction of luminous starbursts and should therefore 
be expected to have a higher fraction of X-ray-detected 
starbursts than would a comparable X-ray selected sam- 
ple. 

2.6.1. Highly Obscured AGNs 

The selection criteria described above will find the bulk 
of the unobscured and moderately obscured AGNs to the 
X-ray survey flux limits. However, the most highly ob- 
scur ed or low-luminosity AG Ns will be mis s ed. F or exam- 
ple, EaddTetin ^QQM) and lFiore et~al1 (|20Q9f ). among 
others, recently found tentative evidence for a population 
of highly obscured AGNs (potentially Ah ~ 10 24 5 cm~ 2 ) 
that lack significant X-ray emission. While the exact con- 
tribution from this highly obscured pop ulation to the to- 
tal number of AGNs is un clear (see e.g.. lAlexander et~al"1 
l2008bt iDonlev et al.ll2008l ). they may represent a numer- 
ically important AGN population. Unfortunately, the 
reliable identification of such AGNs is difficult even in 
nearby sources and is beyond the scope of this work (for 
a comprehen sive analysis and rev iew of infrared selection 
of AGNs, see IDonlev et al.ll2008T) . We instead attempt to 
account for the effect that these missing AGNs have on 
the AGN fraction by adopting th e intrinsic distribu tion 
of column densities determined bv lTozzi et al.l (|2006h and 
comparing it to the observed distribution of the AGNs 
identified in our 70 /jm sample. A detailed discussion of 
our method is given in 21 

3. DERIVATION OF AGN PROPERTIES AND STAR 
FORMATION RATES 

In this section, we describe how we derived various 
AGN properties, such as the intrinsic column density and 
bolometric luminosity, that are required for an accurate 
determination of the AGN fraction. We also describe how 
we estimated the AGN contribution to the observed mid- 
infrared emission, which is required for deriving SFRs for 
these sources. 

3.1. AGN Bolometric Corrections 

The bolometric luminosity is commonly inferred from 
the mean, intrinsic ener gy distribution of a sample of rep- 
resentative AGNs (e.g.. lElvis et al.lll994[ ). The bolomet- 
ric correction, which essentially gives the fraction of the 
bolometric luminosity emerging in a given band, can then 
be determined. Due to the relative ease with which AGN 
X-ray emission may be cleanly measured, the X-ray bolo- 
metr ic correction is commonly used (e.g., Hopkins et al. 
|2007| ). However, a number of complications exist with 
this method. First, knowledge of the intrinsic absorp- 
tion is needed. For bright X-ray sources, the absorbing 
column density can be obtained directly through spec- 
tral fitting. For sources with few counts, the band ra- 
tio can provide a basic estimate of the absorption (e.g., 
lAlexander et al.l l200 5aj) , but some uncertainty will re- 
main. Second, the bolometric correction is known to 
vary with AGN luminosity, due to the luminosity de- 
pendence of the power-law slope between 2500 A and 
2 keV (denoted apx), although th is dependence is well 
quantified fe.g.. iSteffen et al1l2006l) . Third, the intrinsic 
spread in AGN SEDs, possibly due in part to variability, 
is known to result in an uncertainty of a factor of sev- 



eral in the bolometric cor rection (|Hopkins et al.l 120071 ; 
iVasudevan fc Fabian|[200l . 

For the purposes of this study, we use t he luminosity- 
depen dent X-ray bolometric corrections of lHopkins et al.1 
(2007), determined for unobscured (type-1) quasars, 
which transform the intrinsic 2-10 keV luminosity to 
the bolometric luminosity and account for the luminos- 
ity dependence of the SED on «ox, the power-law slope 
between 2500 A and 2 keV. To use this correction, we 
need estimates of the intrinsic luminosity and hence of 
the obscuring column dens i ties fo r the AGNs in our sam- 
ple. Recently, iTozzi et ail (|2006l ) performed X-ray spec- 
tral fitting of the AGNs in the CDF-S to derive reliable 
estimates of the column densit ies. We use t he va lues of 
the column density derived bv ITozzi et al.l (|2006| ) when 
possible. For sources not included in the Tozzi et al. sam- 
ple (80 of 108 AGNs), we estimate the column density 
using the observed X-ray band ratio for each source as 
follows. 

First, we model the X-ray emission using an absorbed 
power-law model in xspec with both intrinsic and Galac- 
tic absorption. In xspec, the model is defined as 
wabsxzwabsxzpow. The photon index of the powe r-law 
component was fixed to 1.8 fe.g-. lTurner et al]|1997l ) and 
the redshift of the zwabs and zpow components was fixed 
to that of the source. We additionally fixed the Galactic 
column density to A H = 8.8 x 10 19 cm~ 2 for the E-CDF-S 
dStark et al II1992D and to Nu = 1.3 x 10 20 cm" 2 for the 
EGS pickev fc Lockma"n1ll990f ). Next, for each source, 
we used this model to find the intrinsic column density 
that reproduces the observed band ratio. To account for 
the Chandra response, we extracted aimpoint responses 
from the event files used to construct the catalogs and 
used the appropriate responses during the fitting. Band 
ratios were taken from the catalogs, except in the EGS, 
where they were calculated from the provided counts and 
adjusted to the aim-point values using the source and 
aim-point effective exposure values. Approximately 10% 
of the X-ray AGNs were not detected in the soft band, 
implying a very hard spectrum and resulting in a lower 
limit to the band ratio. For these sources, we use the 
column density that corresponds to the lower limit. 

Inferred column densities for the AGNs in our infrared- 
selected sample vary from Ah < 10 19 cm~ 2 to Ah ~ 
10 24 cm~ 2 . In genera l, our values agree reasonably well 
with those derived by ITozzi et al . (2006) for a sample of 
194 X-ray-identified AGNs in the CDF-S with detections 
in both the hard and soft bands (a « 1 dex; see Fig- 
ure [4]). Some of the scatter is due to differences between 
the redshifts used by Tozzi et al. and those used by us 
(which have been supplemented with new spectroscopic 
redshifts and the high-quality photometric redshifts dis- 
cussed in 32~3|) . For 53 (« 25%) of the 194 sources, our 
redshifts differ by more than 20% from those used by 
Tozzi et al. Although the agreement between our values 
of the column density and those of Tozzi et al. is good 
overall, at values of Ah > 3 x 10 23 cm~ 2 , our estimates 
of the column density appear to be systematically low. A 
number of these sources were identified by Tozzi et al. as 
Compton-thick. Unfortunately, due to the typically faint 
X-ray fluxes for such sources, the ir reliable identification 
is difficult (e.g.. ITozzi et al.1120061 ) and beyond the scope 
of this paper. 
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Fig. 4. — Comparison of Nu values from Tozzi ct al. (20QH) and 
from our analysis for X- ray-selected AGNs in the CDF-S. The 
source classification from [Tozzi et al . (2006) is indicated by sym- 
bol type, and 70 /im sources are indicated in blue. Sources with 
Nn < 10 20 cm -2 are set to JVjj = 10 2t) cm -2 for plotting purposes. 
A total of 11 sources have this value for both TVjj estimates and 
appear as a single point in the lower left-hand corner. 

Since our sources were chosen to be luminous 70 /im 
emitters, it is possible that the X-ray flux includes a con- 
tribution from star formation. Such a contribution could 
result in derived values of the column density that are 
systematically low. To determine if this effect is likely 
to be significant in the objects in our sample, we fol- 
low the procedure used in i )2.6l to estimate the expected 
contribution from star f ormation using the relations of 
IPersic k Rephaelil ([20071) . As before, we estimated upper 
limits on the star formation rates by assuming that the 
entire observed 70 /im flux is due to star formation. We 
find that the predicted contribution from star formation 
to the observed 2-8 keV luminosity is < 10% in all cases, 
with a mean value of « 0.8%. Therefore, in our AGN- 
selected sample, the AGN emission likely dominates in 
the X-ray band, and emission related to star formation 
is not expected to have a significant effect on the derived 
column densities. This conclusion is supported by Fig- 
ure [4j in which the column densities of AGNs in 70 /im 
sources (shown in blue) do not differ systematically from 
those of the whole population, with the exception of the 
lowest column densities (Ah < 10 21 cm -2 ), where our 
values of the column density do appear to be systemat- 
ically low. However, the corrections required to correct 
the observed fluxes for these low column densities are 
modest and have little effect on our derived bolometric 
luminosities. 

The resulting column densities and the absorbed 
power-law model described above were used to calculate 
corrections to transform the observed-frame 0.5-8 keV 
fluxes to unabsorbed, rest-frame fluxes. To illustrate the 
range of corrections that we find, we plot the corrections 
required to transform from observed 0.5-8 keV fluxes to 
rest- frame 2- 10 keV fluxes as a func tion of redshift in 
Figure [5] (cf. I Alexander et al.l 12008a). The corrections 
generally range from « 0.5-2, but one source with a high 
column density requires a correction of « 5-6. We then 
estimate the bolometric correction required to scale the 
rest-frame 2-10 keV lu minosity to a bolome tric luminos- 
ity from the models of Hopkins et al.l (|2007l ). 

Lastly, we can obtain basic fiducial estimates of the 
black-hole accretion rate from the bolometric luminosity 
by assuming an efficiency for the conversion of the rest 
mass of the accreting material to luminosity. We adopt 
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Fig. 5. — The correction from observed 0.5—8 keV flux to rest- 
frame 2-10 keV flux as a function of redshift for the X-ray-selected 
AGNs. 

an efficiency of e = 0.1, typical of AGNs ac creting at rates 
~ 10% or more of the Eddington rate fe.g. JMarconi et all 
2004). The accretion rate is then M acc = L^ \/(ec 2 ), 
where c is the speed of light. 

3.2. Mid-Infrared AGN SEDs 

Driven by studies of the infrared background, a great 
deal of work has gone into constructing model infrared 
AGN SEDs from the observed SEDs of large samples 
of AGNs. We consider three recent model SEDs to esti- 
mate the AGN contribution to the mi d-infrared emission : 
the type-1 and type-2 AGN models oflSilva et all (|2004j) , 
the type-1 AGN models of lHopkins et al.l (|2007f) . and a 
mea n SEP of AGNs from a flux-limited Sw ift BAT sur- 
vey dTueller et alJl 2008: W inter et all 120091) . 

ISilva et al.l (|2004D constructed type-1 and type-2 AGN 
infrared SEDs using a sample of 33 Seyfert galaxies 
and 11 quasars with available nuclear mid- infrared and 
X-ray fluxes. They constructed intrinsic SEDs by in- 
terpolating the observed SEDs (up to rest-frame A ss 
20 /Ltm). Beyond A w 20 /im, they extrapolated from 
th e observed SEDs u si ng th e radiative-transfer models 
of IGranato fc Dane"sel (|1994D for a number of different 
absorbing column d ensities. In a more recent study, 
iHopkins et all ((20071 ) construct a model SED for type- 
1 quasars using a number of components at different 
wavelengths, including a mean optical spectrum and a 
power-law X-ray spectru m. In the infrared, th ey adopt 
the mean spect r um fro m iRichards et al.l (|2006l ). Lastly, 
iMullanev etaTI (|2009f ) use a sample of 36 AGNs de- 
tected wTth_5£ifeer and th e Swift BAT from the sam- 
ple of iWinter et al.l (|2009|) . selected to have no strong 
indication from PAH features of a significant contribu- 
tion from star formation to the mid-infrared emission. 
These X-ray selected AGNs should be fairly representa- 
tive of the AGNs in our sample. Mullaney et al. have con- 
structed an average mid-infrared SED from these sources 
using a combination of Spitzer IRS spectra and IRAS 
photometry (which provide coverage at wavelengths be- 
yond those covered by the IRS spectra). These three 
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studies provide some of the best determinations of the 
mid-infrared SEDs of AGNs currently available. 

To normalize the AGN SEDs to the observed X-ray 
luminosity, we assume that the mid-infrared SED scales 
linearly with the bolometric AGN luminosity over the 
luminosity range of our sample (i.e., that the AGN lu- 
minosity is the pri nciple determinan t of the infrared lu- 
minosity; see, e.g.. lHaas et al.l 12003). This approxima- 
tio n holds well fo r the lu minosity-dependent AGN SEDs 
of Hopk ins et al.l ([2007) , which include effects such as 
the dependence of aox on the AGN luminosity (e.g., 
ISteffen et al.l [2006). The bolometric luminosity was de- 
rived following q3.ll and scales with the intrinsic X-ray 
luminosity approximately as £^f N oc £2-10 koV over 
the luminosity r ange of our sample. The models of 
iSilva et al.l (|2004[ ) and the average BAT AGN SED that 
we use are luminosity independent (i.e., their shapes do 
not change as a function of AGN luminosity). For con- 
sistency with the Hopkins et al. models, we scaled the 
type-1 AGN Silva et al. models to have the same 2500 A 
luminosity as the (type-1 AGN) Hopkins et al. models at 
a given bolometric luminosity. This same scaling was also 
used when scaling the type-2 AGN Silva et al. models. 
For the BAT AGN SED, the scaling was set such that 
the average Lir of the BAT sample ( calculated from the 
obser ved fluxes using the relation of Sanders fc Mirabel 
1996) is recovered correctly from the bolometric luminos- 
ity corresponding to the average Lx of the BAT sample. 

Figure [6] compares the three AGN models (scaled as 
described above) for AGN bolometric luminosities of 
log(L b oi/£©) = 11-5 and log(L bo i/io) = 12.5. The left 
two p anels show type-1 AG N SEDs and a starburst SED 
(from iCharv fc Elbazi 120011 ) chosen to have roughly the 
same rest-frame 24 /jm luminosity as the AGN SEDs. 
The two AGN SEDs agree to within a factor of 2-3 over 
the infrared region (the Silva et al. models predict higher 
infrared flux out to « 70 /im). As the difference be- 
tween the two type-1 AGN models is smaller than the 
expected systematic errors, we adopt the more recent 
model of Hopkins et al. for subsequent analysis of the 
type-1 AGNs. The starburst SED, while having approxi- 
mately the same rest-frame 24 /im luminosity as the AGN 
SEDs, clearly dominates at longer wavelengths, due to 
the lower temperature of dust the reprocesses emission 
from young stars. 

The right panels of Figure [6] show the av erage BAT 
AGN SED and the type-2 AGN SEDs of ISilva et afl 
(2004) for a variety of Ah values, again with starburst 
SEDs chosen to have approximately the same 24 /im lu- 
minosity as those of the AGNs overlaid. At higher values 
of Ah, the Silva et al. type-2 AGN SEDs show heavy 
extinction at wavelengths below ~ 10 /im compared to 
the type-1 AGN SEDs shown in the left panels. In gen- 
eral, however, between ~ 20-70 /im (the approximate 
range probed by our observed-frame 70 /im data) the 
Silva et al. SEDs are very similar, both among the type- 
2 AGN models of different Ah and when compared to 
the type-1 AGN model. The average BAT AGN SED, 
however, has higher luminosity at wavelengths beyond 
~ 40 /im (by up to a factor of ~ 3 over the probed wave- 
length range) compared to the Silva et al. models. This 
difference will result in larger predicted observed-frame 
70 /im fluxes for type-2 AGNs with z < 0.75 when the 



BAT AGN SED is used. For sources at higher redshift 
(which include most of the high-luminosity sources) , the 
predicted 70 /im fluxes from the two models will agree 
closely. As this difference will primarily affect only the 
lower-luminosity sources (as low-redshift sources tend to 
have lower luminosities), which generally have low pre- 
dicted AGN contributions, our results are not signifi- 
cantly changed from those obtained using the Silva et 
al. models. Therefore, we adopt the Silva et al. models 
for the type-2 AGNs models used in the next section. 

3.3. Predicted AGN contribution to the Mid-Infrared 

Flux 

We can use the model AGN SEDs to predict the AGN 
contribution to the observed 70 /im flux. Before doing 
so, we divided the AGN sample into subsamples based 
on the pre sence of type - 1 AG N optical characteristics in 
either the iBauer et all ()2004D catalog, the COMBO-17 
catalog, or the DEEP2 catalog. In these catalogs, non- 
type- 1 AGNs are simply AGN that are not clearly type-1 
AGNs, so the non- type-1 AGNs are likely to include some 
type-1 AGNs with more subtle characteristics. In total, 
18 (« 17%) of the 108 AGNs in our sample were identified 
as type-1 AGNs. We used type-1 AGN models to predict 
the mid-infrared fluxes of the type-1 AGNs and type-2 
AGN models with the appropriate intrinsic column den- 
sity to predict the mid- infrared fluxes of the non- type-1 
AGNs. We then used the intrinsic, unabsorbed 2-10 keV 
flux derived earlier (see ^3.ip to normalize the models 
(effectively a bolometric correction) and used linear in- 
terpolation (in log space) of the model SEDs to derive 
the observed-frame 70 /im flux or luminosity. In Fig- 
ure [3 we compare the observed 70 /jm luminosity to the 
one predicted by the models. In general, the predicted 
70 /im luminosity is much lower than that observed for 
the majority of AGNs in our sample. 

Critical to this comparison is the normalization of the 
SED, which depends on an accurate estimate of the in- 
trinsic X-ray flux. As discussed in %3.1\ our sample 
of X-ray-detected AGNs does not show evidence from 
X-ray band ratios of being highly extincted (the major- 
ity of inferred column densities are < 2 x 10 23 cm -2 ). 
Therefore, the corrections required to convert observed 
X-ray fluxes to unabsorbed, intrinsic fluxes are typically 
modest and should not be subject to large uncertain- 
ties. Indeed, it is clear from Figure [7] that the predicted 
(observed-frame) AGN 70 /im luminosity (derived from 
the redshifted AGN model) is consistent with the ob- 
served one in all but one system. For this one system (in 
which the predicted luminosity exceeds the observed one 
by a factor of ~ 3), X-ray variability may account for 
the discrepancy, if this AGN was observed in a state of 
higher-than-average X-ray luminosity (conversely, vari- 
ability may also lead to underestimates of the AGN lu- 
minosity in some sources observed in lower-than-average 
states). We emphasize that, because our results depend 
only on properties averaged over many systems, they 
should not be strongly affected by such variability. To 
illustrate this point, in Figure [8] we plot our predicted 
average fractional 70 /im contributions from the AGN 
to the observed flux, binned on 1-dex bins, against the 
observed mid-infrared color for the AGNs in our sample 
with z < 1.5. The mid-infrared color has been found 
to be a rough indicator of the AGN contribution, with 
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Fig, 6. — Comparison o f vario us AGN SED models (denoted by dashed or dashcd-dotte d lines). The typ e-1 AGN models of ISilva et al.l 
(2004) and Hopkins et al. (2007) are shown in the left panels, and type-2 AGN models of ISilva et al. (2004) with ran ge of intrinsic column 
densities are shown in the right panels. Also shown (in red) is the average SED fr om the BAT AGN sample (see £13.21 for details). For 
comparison, the solid green line shows a starburst SED from Chary & Elbaz (2001) with a SFR chosen such that the 24 fim luminosity 
approximately matches that of the AGN SEDs. 



warmer colors indicating higher AGN contributions, at 
least out to z ~ 1.5; beyond this red shift, the color is 
less reliable (e.g. JMurlanev et"aTll2009D . There is a clear 
trend between the two indicators: systems with high pre- 
dicted AGN contributions to the observed 70 /zm flux 
tend to have warmer colors than systems with low pre- 
dicted AGN contributions. Therefore, on average, it ap- 
pears that our method produces estimates of the AGN 
contribution to the mid-infrared flux that are generally 
consistent with dust temperatures indicated by the mid- 
infrared color. Further comparisons to other estimates of 
the relative AGN contribution, such as those from spec- 
tra decomposition of mid-infrared spectra, will provide 
useful tests of systematic errors in our method. How- 
ever, such analyses are beyond the scope of this work. 

Additionally, as the AGNs can contribute significantly 
to the observed mid-infrared emission, sources of a given 
SFR and redshift that host AGNs will be detected more 
readily than those that lack AGNs. To avoid biasing 
our SFR-selected sample toward systems with AGNs, we 
constructed an unbiased sample (henceforth known as 
the "SFR sample") by eliminating sources in which the 
net 70 /jm flux (after subtracting the AGN's contribu- 
tion) results in a signal-to-noise ratio that falls below 
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Fig. 7. — Comparison of the AGN's predicted observed-frame 
70 (im l uminosity to the observed 70 /im luminosity using the 
models of Silva ct al. (2004) for non-type-1 AGNs (circles) and the 
models of Hopkins ct al. (2007) for type-1 AGNs (triangles) . Open 
symbols denote systems that have a net S /N < 3 after subtracting 
the estimated contribution from the AGN (see text for details). 
The solid line denotes equality, and the dashed and dotted lines 
indicate an AGN contribution to the observed 70 fim luminosity of 
1% and 10%, respectively. 

our adopted limit (S/N = 3; see g2j) . We found that 23 
(20%) of 108 AGNs fell below this limit in the combined 
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Fig. 8. — The observed mid-infrared color versus the ratio of 
(observed-frame) predicted-to-observed 70 ^m flux for the AGNs 
in the sample with z < 1.5. The filled points show the mean val- 
ues (calculated using the non-logarithmic values of the colors) of 
subsamples of objects in bins with widths indicated by the hori- 
zontal error bars. The open points show the median values for the 
same bins. Vertical error bars indicate the standard deviation of 
the colors in each bin. 

E-CDF-S and EGS samples (shown as open symbols in 
Figure [7]) . The SFR sample is used only to study the 
AGN activity as a function of SFR. When we examine 
the AGN activity as a function of other properties (e.g., 
mid- infrared color), the full sample of 108 AGNs is used. 
We note that most of the AGNs that were eliminated 
have log (F24 IJ , nl /F- 70 ^ m ) > —0.7 and are predicted to con- 
tribute a large fraction (> 50%) of the 70 /mi flux. In 
the SFR sample of 85 AGNs (of 1022 70 jum sources in 
total), only two AGNs have an estimated contribution 
to the observed 70 /im flux of more than 50%. There- 
fore, this cut eliminates most of the sources for which 
the determination of the SFR is likely to be subject to 
large systematic uncertainties (i.e., those in which AGN- 
powered emission likely dominates at 70 /im). 

Lastly, due to intrinsic differences of the SEDs of 
systems with the same SFR, some systems at a given 
SFR will have 70 /im fluxes that fall below our adopted 
flux limit, resulting in our sample being incomplete 
at the given SFR. We correct for this incompleteness 
by estimating t he sc atter about the relation used by 
iCharv fe Elbazl ()200lD for local starbursts. We estimated 
a scatter of ~ 0.5 dex in the observed flux at 30-70 /im 
at all star formation rates, and we have used this value 
to estimate the likely number of missed starbursts for 
each detected one. We did this by generating, for each 
detected source, a normal distribution of fluxes around 
the observed 70um flux. We then calculate the fraction 
of sources that fall below the flux limit for the detected 
source under consideration and correct the total num- 
ber of starbursts by the sum over all sources. We as- 
sume that these missed sources, since they are presum- 
ably cooler than the average starburst for a given LIR, 
do not host AGNs. Again, as in the rest of our analysis, 
we do not include any evolution in the starburst SEDs, 
as such evolution is currently poorly understood. This 
correction is generally small, and is significant only for 
those sources detected near the flux limit. Additionally, 
where the AGN fraction is high (such as at high SFRs), 
the effect of this correction is reduced (since the frac- 
tion of non-AGN hosts, which is being adjusted, is by 
definition lower). 



3.4. Star Formation Rates 

As discussed in §1, mid-to- far- infrared observations 
sample a significant fraction of the energy emitted by 
massive stars in dusty environments, with much of the 
remainder emerging at UV wavelengths. Therefore the 
total star formation rate may be estimated as the sum of 
the rate inferred from direct UV emission and the rate 
inferred from the reprocessed, infrared emission (e.g., 
IDaddi et al.l[2f)07ah . 

To trace the infrared emission associated with star 
formation, we use the observed-frame 70 /mi luminos- 
ity. A number of recent studies (e.g. , iShi et all l2007t 
iTadhunter et al.ll2007t IVega et al.ll2007j ) have found that, 
at a given star formation rate, the observed-frame 70 /mi 
emission suffers from significantly less AGN contamina- 
tion and spectral complexity than 24 /tm or shorter- 
wavelength emission, particularly at higher redshifts (z > 
1.5). However, the contribution from AGNs to the 70 /tm 
emission can still be significant. Therefore, for sources 
hosting an AGN, we estimated the AGN contribution to 
the 70 /mi luminosity using empirical AGN SEDs (see [J3] 
for details). The net observed-frame 70 /mi luminosity 
was then converted to a rest-frame, 8-1000 /im luminos- 
ity (denoted £ir) usi ng the dusty starburst models of 
ICharv fc Elbad (|2001[ ). which are luminosity depe ndent, 
and the prescription of iSanders fe Mirabell ((T996). The 
resulting Z /tr was th e n con verted to a SFR using the 
relatio n of iKennicuttl (JT998) , which assumes a iSalpeterl 
(fl955l) IMF, as follows: 



SFR 



IR 



M Q yr- 



1.73 x 10 



-10 



We also investigated the use of other publicly available 
dusty starburst models and found that, for our sample 
and using the observed-frame 70 /mi flux (which samples 
rest- frame wavelengths > 20 /tm), there is little practical 
difference (< 50% in -Lir) betwee n the Chary fe Elbaz 
models that we adopt and those of IDale fe Heloul (2002) 
or lRieke et all (|2009D Fl We note that all of these models 
are derived from local samples of starburst galaxies and 
hence could differ systematically from the sources in our 
stud y (which are genera lly at redshifts of 0.5-1.5). How- 
ever, lElbaz et all (12002) in a study of infrared-luminous 
galaxies, found good agreement out to z ~ 1 between the 
radio-derived SFRs and those derived f rom the infrared 
emission. At z ~ 2. IDaddi et all (j2007bl ) also find reason- 
able agreement between various indicators of star forma- 
tion rate, including the infrared (we note, however, that 
sub-millimeter galaxies, at z > 2, likely have lower dust 
temperatures for a given SFR than local ULIRGs, e.g., 
iPope et al.ll2"0"07t iCoppin et ai1 l2008) . Since the choice of 
model makes little difference for the values of Lir that 
we derive and a clear consensus as to the most appropri- 
ate model for high-redshift, high-luminosity sources has 
yet to emerge, we adopt the models of Chary & Elbaz 
for all further analysis. 

As stated above, an additional important component 
of the bolometric emission from star formation emerges 

16 If, however, we were to use the observed-frame 24 /im flux, 
which would sample rest-frame wavelengths of > 5 fira for our 
sample, the differences between models can become large (up to 
factors of > 10 in Ltr), as illustrated by Salim et al. (2009). 
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in the UV. This "direct" emission must be included 
when deriving the total SFR. To this end, the rest-frame 
UV luminosity was estimated using the rest-frame B-, 
V-, and i?-band fluxes derived from the observed SED 
by fitting galaxy templates using ZEBRA or by simple 
linear interpolation (when pos sible). The UV co nver- 
sions from two r ecent studies, Daddi et al.l f|2004T ) and 
IBell et al.1 d200l . were used to transform the UV lu- 
minosity to a SFR. The Daddi et al. relation uses the 
1500 A rest-frame luminosity to calculate the UV SFR 
as SFR 15 oo/(M vr" 1 ) = 1.13 x 10- 28 (L 150 o/L®). Bell 
et al. use the 2800 A luminosity to estimate the UV 
SFR as SFR 28OO /(M yr" 1 ) = 8.99 x 10 - 29 (i 2800 /L Q ). 
At z < 1.5, the 1500 A rest-frame emission is sam- 
pled only by GALEX observations. Therefore, to avoid 
large extrapolations for sources at z < 1.5, we use the 
Daddi et al. relation only for those sources that have 
GALEX detections (« 40%). The Daddi et al. and 
Bell et al. estimates, which typically agree to within a 
factor of a few, were then averaged to obtain the UV 
SFR, SFR UV . The total SFR is then calculated as 
SFR to t = SFR m + SFR UV . It should be noted that 
no correction is applied to account for extinction in the 
UV, as emission that is absorbed by dust will be repro- 
cessed and is therefore included in the infrared-derived 
SFR. Since emission from the AGN may dominate at UV 
wavelengths, we assumed that the fraction of emission 
emerging in the UV from star formation for the AGN 
sources is the same as that for the non-AGN sample on 
average (i.e., the extinction in the UV is similar). For 
the sample as a whole, the UV SFRs generally represent 
< 50% of the total SFR. Figure [2]: shows the distribution 
of SFRs. It is clear from this figure that the sample is 
approximately complete to z — 1.0 (z = 2.0) for sources 
with SFR > 100 Mq yr" 1 (SFR > 600 M© yr" 1 ). 

4. CALCULATION OF THE AGN FRACTION 

The AGN fraction is defined as the number of AGNs 
above a given intrinsic X-ray luminosity divided by the 
total number of sources in which an AGN was detected 
or could have been detected, g iven the sensitivity limits 
of the X-ray observations (e.g.. lLehmer et af] |2008) . The 
cu mulative AGN fra c tion m ay then be calculated follow- 
ing iSilverman et al.l (|2008f ) so that the contribution of 
each AGN to the total fraction is included: 



Wagn 



/agn - ^2 



1 



i=l 



A 



al , i 



(4) 



In this equation, Aagn is the total number of AGNs in 
the sample and A ga i,i is the number of galaxies in which 
the i-ih AGN could have been detected. We further re- 
strict Nggxj to include only those sources that lie in re- 
gions of sensitivity great enough to detect (at a S/N > 3) 
the 70 fim flux of the z-tli AGN, thereby imposing a flux 
limit (as opposed to a S/N limit) on the sources that 
contribute to each AGN's contribution to the total frac- 
tion. Th e error in the AGN fra ction is calculated (again 
following Silv erman et al.l [200 8 ) as: 



(5) 




where cr p hy S is the contribution to the error from uncer- 
tainties in the physical properties used to define the bins 
(the SFR, the rest-frame mid-infrared luminosity, and 
the mid-infrared color) and o ^v H is the uncertainty re- 
sulting from the probabilistic treatment of the intrinsic 
Ah distribution (discussed in detail later in this section). 
The (Jphys term is estimated using a Monte Carlo tech- 
nique as follows. For each source, we drew random val- 
ues of the physical property from a normal distribution 
centered on the measured value of the physical property 
with a standard deviation given by the uncertainty in 
the property (e.g., for the SFR, we used errors derived 
from the reported errors in the 70 /im fluxes). We re- 
peated this procedure 100 times, each time calculating 
a new fraction, and estimated a p hy S from the resulting 
distribution of /agn • 

Because the sensitivity of the X-ray data used in this 
study varies with position (by factors of > 10), system- 
atic errors will be induced in the AGN fraction if this 
variation is not accounted for when determining A ga i,.;. 
To remove the effects of X-ray sensitivity variations, we 
include in Ngaij only those galaxies in which an AGN 
with luminosity Lx,i could have been detected if present 
(i.e., only galaxies with limiting X-ray luminosities below 
Lx.i)- To estimate the X-ray sensitivity limits, we used 
sensitivity maps generated separately for each field. Due 
to the dependence of the Chandra point spread function 
and effective area on the off-axis angle, the X-ray sen- 
sitivity across a single field is a strong function of the 
position relative to the aim point of the observations. 
This variation can be estimated and, under the a ssump- 
tion of Poisson statistics (e.g., iLuo et all |2008|) . maps 
may be generated that give the sensitivity limit of the 
survey as a function of position. For the CDF-S and 
E-CDF-S, maps were generated as part of the cata log 
construction (see lLehmer et aLll2005HLuo et al.ll2008t l in 
terms of the limiting flux that corresponds to the num- 
ber of counts required for the secure detection of a source 
with a r = 1.4 power-law spectrum. For the EGS, maps 
were provided directly in terms of limiting counts for a 
source with the same spectrum (see Laird et al. 2009). 
If a source lies in the CDF-S region (and therefore has 
both 2 Ms and 250 ks coverage), we adopted the lowest 
limiting flux from the CDF-S and E-CDF-S sensitivity 
maps at that position. This flux is generally that of the 
2 Ms CDF-S except in some regions at large (> 8') an- 
gles from the average CDF-S aim point that lie near the 
E-CDF-S aim points. 

Next, we attempt to account for the distribution of 
intrinsic AGN column densities, which will affect both 
the overall AGN fraction (as high-column-density sources 
will tend to be missed in even the deepest X-ray sur- 
veys) and will result in field-to-field variations in the 
fraction due to different exposure times that probe dif- 
ferent column densities at a given redshift. For example, 
at a given luminosity and redshift, AGNs with larger 
intrinsic column densities can be detected near the cen- 
ter of the CDF-S (where exposure times reach rj 2 Ms) 
than in the shallower E-CDF-S (where typical exposure 
times are « 250 ks). We account for these effects by 
using an estimate of the intrinsic distribution of AGN 
column densities to effectively adjust the detection lim- 
its derived from the sensitivity maps. We use the in- 
trinsic AGN Ah distribution determined for the CDF-S 
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by iTozzi et al.l (|2006[ ) , who found that the distribution 
can be approximated by a log-normal distribution with 
<log^ H /(cm- 2 )) w 23.1 and er « 1.1. 

In deriving this distribution, Tozzi et al. assumed there 
is no strong dependence on the distribution with intrinsic 
luminosity or redshift. However, there is some evidence 
that the absorb ed fraction is lower for higher luminos- 
ity AGNs (e.g. , lUeda et al.ll2Q0l iTreister fe UrTvl [20051: 
lHasingerll2008D . Given the current uncertainties in the 
detailed dependence of the Ah distribution on luminosity 
and redshift, we do not attempt to account for any such 
dependence but note that the Tozzi et al. distribution 
was derived using an X-ray-selected AGN sample that 
is similar to our AGN sample (and shares our CDF-S 
sources) and should therefore apply well to our sample 
on average. We slightly modified the log-normal distri- 
bution d escribed above to match better the actual one 
found bv lTozzi et al.1 ((2006) by maintaining a flat distri- 
bution from A^h ~ 10 23 cm" 2 to Ah ~ 10 24 cm" 2 and 
by including the 10% of objects at low values of .Ah 
(< 10 20 cm" 2 ). We note that this distribution is also 
generally co nsistent with that adopted in other re cent 
studie s (e.g..lGilli et al.ll2007l:lMerloni fc Heinzll2008f ) and 
in the Hopkin s et al.l (|2007f ) study from which we have 
derived the AGN bolometric corrections. 

Using the sensitivity maps described above, we can 
now determine the number of 70 /im sources in which 
an AGN of a given luminosity, subject to the adopted 
A^h distribution, could have been detected. We first 
draw 1022 values of Ah from the above distribution 
and assign these values to each source. Then, for the 
ith AGN in our sample, we place hypothetical AGNs 
with intrinsic 0.5-8 keV luminosities equal to Lx,i in- 
all 70 /im sources. We then calculate the resulting ab- 
sorbed, observed-frame 0.5-8 keV fluxes for each hypo- 
thetical AGN using the source redshifts, the assigned col- 
umn densities, and the absorbed power-law AGN model 
described in §3.11 Although this model does not in- 
clude the scattered, reflected, or lin e emission identified 
in many highly obscured AGNs (e.g.. lMalizia et aLll2003l) 
and hence may underpredict the soft flux emerging from 
sources assigned high values of A^h (above ~ 10 24 cm -2 ), 
we found no significant difference in our results when 
an em pirically motivated AGN model (constructed fol- 
lowing lAlexander et al.l[2005b| ) was used. Therefore, for 
simplicity, we adopt the absorbed power-law model in 
our analysis. 

We can now compare the predicted observed flux for 
each hypothetical AGN with the sensitivity limit at its 
position to determine whether or not such an AGN could 
have been detected if present. Since the Chandra re- 
sponse is a strong function of energy, and AGNs with 
higher values of A^h have harder spectra, the resulting 
number of predicted 0.5-8 keV Chandra counts (which 
determines whether or not a source will be detected) will 
differ for two sources with similar observed-frame 0.5- 
8 keV fluxes and redshifts but different values of A^h- 
Therefore, for each 70 /im source, we calculate the num- 
ber of detected counts expected from a given hypotheti- 
cal AGN using the absorbed power-law AGN model (at 
the appropriate redshift), the appropriate Chandra re- 
sponse, and the effective exposure time at the position 
of the source. We then compare the predicted number 



of counts to the limiting number of counts at the source 
position. If the predicted number of counts exceeds the 
limiting number, the source is included in the calculation 
of the AGN fraction. We use the full-band (0.5-8 keV 
for the CDF-S and E-CDF-S and 0.5-7 keV for the EGS) 
limiting counts given by the sensitivity maps described 
above. For the CDF-S and E-CDF-S, for which the maps 
are given in terms of a limiting flux calculated assuming 
a r = 1.4 power-law spectrum, we convert the flux to 
counts using the same spectrum and the effective expo- 
sure times at the positions of the sources. Finally, we re- 
peat this entire procedure 100 times, and adopt the mean 
value of the fraction, < /agn > , as the best estimate and 
set <7at h to the standard deviation of the fraction over the 
100 runs. 

To illustrate the effect of including the column den- 
sity distribution in the calculation of the AGN fraction, 
we plot in Figure [9] the cumulative AGN fraction as 
a function of the AGN luminosity for the three fields. 
Our method of accounting for the column density dis- 
tribution increases the AGN fraction overall, but par- 
ticularly at lower redshifts and hence, on average, at 
lower AGN luminosities (by « 10-30% at L . 5-8.0 keV ~ 
2 x 10 41 erg s" 1 ). At the highest redshifts and AGN lumi- 
nosities in our sample (L . 5-8.0 keV ^ 10 44 erg s" 1 ), the 
fraction is altered only slightly, since the emission probed 
by the Chandra data is at higher rest-frame energies and 
hence less affected by the obscuration. 

Additionally, the correction produces larger changes 
in the shallower fields, such that the AGN fraction in 
the CDF-S, the field with the deepest X-ray data, has 
the least change and the EGS the greatest. This result 
is expected given that a larger fraction of high-column- 
density sources will be missed in the shallower X-ray 
fields at a given intrinsic AGN luminosity (and hence 
the AGN fraction will be biased low). We note, how- 
ever, that, even after the correction is applied, the EGS 
tends to have the lowest cumulative AGN fraction at a 
given AGN luminosity (particularly at lower luminosi- 
ties), suggesting that our adopted column-density distri- 
bution may not apply as well to the sources in the EGS 
as to those in the E-CDF-S (possibly due to cosmic vari- 
ance) or that there may be some evolution in the cumula- 
tive AGN fraction with redshift (as the EGS data probe 
lower redshifts on average than the E-CDF-S data at a 
given AGN luminosity). Despite this issue, the cumu- 
lative fractions of the three fields are roughly consistent 
with one another given the uncertainties. Therefore, to 
obtain a larger sample size, we henceforth examine the 
AGN fraction of the combined E-CDF-S and EGS sam- 
ple. 

5. RESULTS 

As discussed in §1, the AGN fraction, which gives the 
detection rate of AGNs in a given sample, is related to 
the duty cycle of AGN activity. Higher fractions im- 
ply that the AGNs in these systems spend more time in 
active states than do AGNs in systems with lower frac- 
tions. Therefore, the AGN fraction gives an indicator 
of the ubiquity of black-hole growth. Along with esti- 
mates of the relative levels of bulge and SMBH growth, 
this information can be used to understand how SMBHs 
grow during periods of vigorous starburst activity. In 
this section, we investigate the relation between SMBH 



14 



Rafferty et al. 



1 .000 E 



0.100 



0.010 r 



0.001 




42 43 44 

l°gL . 5 -8kev (ergs -1 ) 



45 



1.3 

i 1.2 

z" 
o 
< 

" 1.1 

z 
o 
< 

— 

1.0 

0.9 
41 




42 43 44 

l°g L o. 5 -8kev (ergs" 1 ) 



45 



Fig. 9. — Left: The AGN fraction corrected for the Tozzi et al. (2006) distribution of column densities (/agn) as a function of the 
intrinsic, rest-frame 0.5—8 keV AGN luminosity. Right: The ratio of the corrected-to- uncorrected AGN fraction as a function of the AGN 
luminosity. The fraction is plotted separately for each field, as indicated by the different symbols (note, however, that the E-CDF-S sources 
include those of the CDF-S). Errors are calculated following equation (0, but account only for the sampling error (i.e., <r p hys = and 
°-JVh = 0). 



growth and the SFR for our sample of starbursts, and 
we examine the dependence of the AGN fraction on the 
observed source properties. Properties of interest include 
the mid-infrared color (which gives a rough measure of 
the temperature of the emitting dust and thus the rela- 
tive contributions from AGNs and star formation to the 
power source of the reprocessed emission), the rest- frame 
mid-infrared luminosity, and the star formation rate. 

For each of these properties, we use the com- 
bined E-CDF-S and EGS 70 /im sample and inves- 
tigate two minimum rest-frame AGN cutoff luminosi- 
ties: L . 5-8.0 kov > 10 41 erg s" 1 and Lo.5-8.0 koV > 
10 43 erg s _1 . However, due to the flux-limited nature of 
the X-ray and far-infrared surveys upon which our anal- 
ysis is based, both this cutoff luminosity and the rest- 
frame mid-infrared luminosity and SFR are increasing 
functions of redshift (see Figured]). Therefore, caution 
must be exercised when interpreting trends in the AGN 
fraction with luminosity or SFR when a large range in 
cutoff AGN luminosities is present (for reference, the cut- 
off AGN luminosities are indicated on the relevant plots). 

When examining the AGN fraction as a function of 
SFR, we used the subsample of 85 AGNs (and 999 70 //m 
sources in total) created by filtering out sources that fall 
below our adopted 70 /im S/N limit after subtracting 
the estimated contribution from the AGN to the 70 /im 
flux (see ^3.3p . To avoid situations in which a single 
AGN dominates the fraction in a given bin, we con- 
struct the bins so that each contains a minimum of 10 
AGNs and exclude any bin in which a single AGN con- 
tributes 30% or more to the fraction in that bin (the 
contribution of each AGN to the total fraction depends 
on the distribution of limiting luminosities; see equation 
[4])- This method minimizes the effect on the fraction 
of a single AGN that might have, for example, an in- 
correct redshift estimate (such sources are expected to 
account for < 10% of our sample; see §2.3j) . Addition- 
ally, the AGN fraction is strictly valid only when cal- 
culated for complete samples (e.g., for all galaxies with 
10 < SFR < 30 Mq yr- 1 and redshifts less than the 
limiting redshift for a SFR = 10 M@ yr -1 galaxy). For 
our adopted minimum number of AGNs per bin (10), 
we found that the typical bin width is small enough such 



that the difference in limiting redshifts at the bin bound- 
aries is much smaller than the typical limiting redshift of 
objects in that bin. Therefore, the samples in each bin 
should be roughly complete. 

5.1. The AGN Fraction and Mid-Infrared Color 

We begin by showing in Figure [TU] the AGN fraction as 
a function of the mid-infrared color of our 70 /im sources. 
It is clear from this figure that the fraction is a strong 
function of the mid-infrared color, rising from 5-10% at 
the smallest values of the ratio (-F24/-F70) (indicative of 
cooler dust temperatures) to ~ 60-70% at the largest 
values. The traditional dividing point of this ratio for 
z ~ sources between dust powered by AGN-dominated 
emission and that powered by star forma tion-dominated 
emiss ion is at logfJas/JR o) ~ —0.7 (e.g.. Ide Grijp et all 
[19851: ISanders et al.Nl988l) . This division also occurs at 
roughly the same ratio when Spitzer bands are used (i.e., 
log[i<24/-FYo] —0.7) and appears to hold out to at least 
z ~ 1.5 (e.g- JMullanev et al .1120091) . Indeed, at approx- 
imately this ratio, the AGN fraction appears to reach 
its maximum, implying that ~ 60-70% of such sources 
host an AGN. Below log(F 2 4/F 7 o) « -0.7, the fraction of 
sources of a given color that hosts an AGN falls rapidly as 
the color indicates cooler temperatures. Although care 
must be taken in interpreting the color at z > 1.5 due 
to the increasing complexity of typical AGN and star- 
burst spectra at the rest-frame wavelengths sampled by 
the 24 /im band in particular, th is result extends the 
analyses of iMullanev et ail (|2009() by showing that the 
mid-infrared color is a useful indicator of luminous AGN 
activity in the distant universe. 

5.2. The AGN Fraction and Mid-Infrared Luminosity 

In Figure [TTJ we plot the AGN fraction against the 
rest-frame 30 /im luminosity of the source (note that 
this quantity includes any AGN contribution) . It is clear 
from this figure that the fraction of sources hosting an 
AGN depends on the rest-frame mid-infrared luminos- 
ity, with higher fractions in sources with higher lumi- 
nosities. The dependence becomes stronger when lower- 
luminosity AGNs are excluded, such that the fraction of 
sources hosting an AGN with Lq. 5-8.0 kcV > 10 43 erg s _1 
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the mean width of which is indicated in the lower right-hand corner. Variations in the fraction on scales smaller than this bin width are 
not significant. The shaded region indicates the 1-cr errors. 



rises with the mid-infrared luminosity from a few per- 
cent at L 30 « 10 10 L Q to - 60% at L 3Q w 5 x 10 12 L Q . 
Therefore, more luminous mid-infrared sources, such as 
ULIRGs, are ~ 10 times more likely to host a luminous 
AGN than lower-luminosity sources, such as local star- 
bursts. Such a result is to be expected if the AGN con- 
tributes to the mid-infrared emission, and the average 
contribution increases with increasing AGN luminosity 
(as is suggested by Figure [7]). 

The peak in the AGN fraction (/ AGN - 50-60%) in 
Figure [TT] at the highest mid-infrared luminosities (cor- 
responding broadly to ULIRG luminosities) is consistent 
with the fraction fa 60%) of loca l ULIRG s at these 
SFRs identified by lAlexander et al.l (|2008aD as hosting 
X-ray AGNs or classified as Seyfert galaxies by iVeilleuxl 
(2006). It is also consistent with the total AGN frac- 
tion of ~ 50-65% implied by the detection rate of highly 
obscured AGNs (identified using Spitze r spectroscopy) 
in a s ample of local ULIRGs studied by llmanishi et al.l 
(|2007l) . We note, however, that the luminosity cutoffs in- 
dicated in Figure QT] do not strictly correspond to LIRG 
and ULIRG cutoffs since, for example, cooler ULIRGs 
may have rest-frame 30 /im luminosities below the indi- 
cated cutoff. Furthermore, the differing selection criteria 
of these studies makes it likely that they probe AGNs 
of different luminosities on average (as noted earlier and 
indicated in Figure HT1 the minimum AGN luminosity to 
which we are sensitive varies with the redshift and hence 
rest-frame 30 jim luminosity). Howe ver, if we restrict 
the sample of lAlexander et all ()2008al ) to those systems 
hosting an X-ray AGN with L2-10 kcV > 10 43 erg s" 1 , 
which matches our selection criteria more closely, we find 
that 4 of 10 local ULIRGs host such an AGN, with an 
implied AGN fraction of ~ 40%, in good agreement with 
our result. Therefore, our results suggest that the AGN 
fraction in the distant ULIRGs in our sample is similar 
to that of local ULIRGs and, furthermore, that the rela- 
tion between AGN activity and star formation in these 
two populations is similar. 

5.3. The AGN Fraction and Star Formation Rate 

Figure [T^] shows 
for AGNs with L .5 
£0.5-8.0 kcV > 10 43 



the cumulative AGN 
> 10 41 erg s 



fraction 
1 and 



Before calculating a SFR, we subtracted the expected 
AGN contribution from the observed 70 /im flux as de- 
scribed in i )3.3l We then used the procedure described 
in M3-4I to infer a SFR from a given net 70 fim flux. 
It is clear from this figure that the AGN fraction in 
sources with lower SFRs (SFR ~ 10 M Q yr _1 ) de- 
pends strongly on the AGN luminosity (£0.5-8.0 kev)- 
When lower-luminosity AGNs are included, the fraction 
of sources hosting an AGN at low SFRs is quite high at 
~ 10-20%. When only higher-luminosity AGNs are con- 
sidered (£0.5-8.0 kcV > 10 43 erg s _1 ), the AGN fraction 
at low SFRs is much lower, rising from w 4-5% at SFRs 
of roughly 10 M@ yr^ 1 to ~ 40% at the highest SFRs as 
/agn oc SFR 75 . 

Additionally, SMGs, which generally lie at high red- 
shift (z ~ 2) where the reprocessed emission from 
cool dust peaks at observed-frame submillimctcr wave- 
lengths, have larger average SFRs than most of the 
sources present in our sample. Their typical SFRs 
are on the order of 1000-3000 M Q yr _1 , and the 
reprocessed emission in these system s appears to be 
largely powered by star formation ( e.g. .lAlexander et al ' 
2005al iMenendez-Delmestre et al.l 120071 iValiante et al 



8.0 kcV 

erg s _1 as a function of SFR 



I2007t iPope et all 12008). SMGs also appear to have a 
high AGN fractio n (~ 20-50%. lAlexander et al] 120051* 
lLaird et"aT1l2010| ). and are thought to be the short-lived 
phase of rapid accretion during a massive merger. Their 
relation to local ULIRGs is unclear (e.g., iSaiina et all 
l2007t iSvmeonidis et all I2009D . but the two classes ap- 
pear to share many of the same properties, with SMGs 
being roughly scaled up by an or der of magnitude 
in far- infrared lumin osity and SFR (Sa una et all 120071 
iVahante et alll2007l ). We note that although our sam- 
ple does not contain a sufficient number of SFR > 
1000 M Q yr _1 objects to determine a reliable AGN 
fraction in this regime, extrapolation of our determi- 
nation of the AGN fraction (Figure [12} to SFRs of 
~ 2000 M Q yr _1 results in an AGN fraction roughly 
consiste nt with the frac t ions f ound by lAlexander et all 
(|2005bf ) and lLaird etafl (j2010Q for SMGs using indepen- 
dent methods and a different field (the CDF-N). 

To investigate whether the AGN fraction has a depen- 
dence on redshift, we plot in Figure [TBI the AGN fraction 
against the SFR for two subsamples divided by the me- 
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dian redshift of the AGN sample, z = 0.9. Due to the 
sensitivity limits of the FIDEL survey, at high redshifts 
only the most luminous 70 /jm sources (with correspond- 
ingly high SFRs) are detected. Conversely, at low red- 
shifts, no very luminous 70 //m sources are present in 
the survey fields. Due to this effect, the two subsamplcs 
do not have sufficient overlap to judge whether there is 
any dependence of the AGN fraction with redshift at a 
given SFR. Additionally, as discussed earlier, the mini- 
mum AGN luminosity for which the cumulative fraction 
is calculated changes with the SFR, as indicated in Fig- 
ure 1131 This changing minimum is largely responsible 
for the trend of decreasing AGN fraction with increas- 
ing SFR seen at z < 0.9 in the left panel of Figure [T51 
Two opposing effects sum to create this trend: (1) the 
number of lower-luminosity AGNs per bin decreases with 
increasing SFR due to a steady rise in the minimum cut- 
off AGN luminosity with increasing redshift (and hence 
SFR) and (2) the contribution to the cumulative fraction 
from higher-luminosity AGNs increases with the SFR (a 
trend visible in the right panel). Overall, the net effect is 
a gradual decrease in the AGN fraction with increasing 



SFR at z < 0.9. At z > 0.9, the change in the cut- 
off AGN luminosity across the sampled range of SFR is 
small, and should not have a large effect on the overall 
trend. 

To examine the effects that AGNs with warm mid- 
infrared colors (i.e., those sources that may have large 
AGN contributions to their 70 fim fluxes) have on our 
results, we have further filtered our sample to exclude all 
sources with \og(F2i / F^) > —0.7. The only significant 
effect of this filtering is reduce the AGN fraction overall 
(and to reduce the sample such that a smaller range of 
SFRs is probed) . The same trends are present both with 
and without this filtering. 

5.4. Black-Hole and Bulge Growth 

It is now well established that black holes and 

their host bulg es are intimately connected (e.g., 

iFerrarese fc Fordl 12005). To illustrate how the growth 
rates of the AGNs in our sample compare to those of 
their bulges, we binned our AGN sample by SFR and 
calculated the mean bolometric AGN luminosity (calcu- 
lated following i )3.ip in each bin. In this comparison, we 
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implicitly assume that the bulge growth rate is approxi- 
mated by the total SFR, although this may not be strictly 
true (e.g., in late-type galaxies). In the left panel of Fig- 
ure [TH we compare these luminosities to the 8-1000 /im 
luminosities from star formation derived following equa- 
tion ©. 

It appears that, in systems identified as having AGNs, 
the SMBHs and bulges in our sample are growing concur- 
rently on average, across a wide range of SFR (roughly 
two orders of magnitude), and at relative rates that 
would produce or maintain the scaling observed lo- 
cally (calculated using the conversions described above 
and assuming SFR /A/bh oc Mbuigc/MBH 700; e.g., 
lHaring fe"~ Rix 2004). We note that we are likely missing 
many systems with low SFRs that would fall mostly in 
the lowest bin. Depending on the average AGN lumi- 
nosity of such systems, our value for the average AGN 
luminosity in this bin could be biased either high (if 
the missed systems have weak AGN activity on average) 
or low (if they have strong AGN activity on average). 
The effect of systems in the latter category (with high 
-f'AGN, boi and low Lsf, ir) will be reduced somewhat 
due to the global increase in specific SFR with redshift, 
which will tend to reduce the fraction of luminous AGNs 
that lack significant star formation. Nevertheless, the un- 
certainty in the lowest bin is large and its value should 
be treated with caution. 

However, because many galaxies at a given SFR lack 
luminous AGN activity, the average AGN luminosity 
(and hence SMBH growth rate) across all galaxies in a 
given bin will be somewhat lower. We can estimate this 
average AGN luminosity by multiplying through by the 
appropriate AGN duty cycle (traced by the AGN frac- 
tion shown in Figure fT2|). This estimate is shown by the 
solid line in the left panel of Figure Q3] and is well below 
that expected from the local scaling relation (with the 
exception of the lowest-SFR bin, which, as noted above, 
is likely strongly affected by incompleteness), implying 
that the bulges in our sample could be growing at a faster 
rate on average relative to their SMBHs than expected 
from the local relation. However, we note that the bulge 
growth rate will often be significantly less than the total 
SFR of the galaxy and, for this reason (among others dis- 
cussed in the next section), this average AGN luminosity 
should be interpreted as an approximate lower limit to 
the true one. 



Additionally, in the right panel of Figure HH we show 
the distribution of the ratio of the SFRs of the host 
galaxies to the SMBH accretion rates calculated in ij3.ll 
Again, it is clear that the ratio of SFR to accretion rate 
for our sample of AGNs is broadly consistent with that 
expected from the median observed ratio of bulge mass 
to black-hole mass found locally, although our sources 
tend to lie at slightly higher ratios (median SFR/Mbh ~ 
1000). However, given the uncertainties in the calcula- 
tion of these quantities, our results agree well with ex- 
pectations from the MBH-A^buige relation and studies of 
the volume-w eighted bulge an d SMBH growth in local 
galaxies (e.g... iHeckman et"aT1l2004l ). 

In both panels of FigurefMl there is large scatter about 
the local scaling. Much of this scatter may be due to the 
systematics discussed above, but could also indicate that 
some of the systems are undergoing large deviations (by 
factors of up to ~ 100) from concurrent SMBH-bulgc 
growth. 

6. DISCUSSION 

In the previous section, we presented the results of 
an analysis of SMBH growth in distant luminous star- 
bursts. We investigated the AGN fraction as a function 
of a number of physical properties and examined the rel- 
ative growth rates of the SMBHs and their host galaxies. 
In particular, we found a strong dependence of the AGN 
fraction on the SFR in our sample, with the fraction ris- 
ing to ~ 30% at the highest SFRs, implying a high duty 
cycle of luminous AGN activity in such systems. 

However, as discussed previously, due to the sensitiv- 
ity limits of the X-ray data and the fact that high-SFR 
objects tend to lie at higher redshifts, we are not gener- 
ally sensitive to lower-luminosity AGNs in these systems. 
Therefore, the total cumulative fraction of AGNs in in- 
termediate and high-SFR systems may be higher than 
shown in Figure [121 Additionally, as discussed in §2.6-H 
it is likely that we have missed a significant population 
of AGNs with very high column densities due to their 
weak X-ray emission. Although we have attempted to 
correct for this missing population using an estimate of 
the distribution of intrinsic AGN colu mn densities, our 
adopted distribution (based on that of iTozzi et al.ll2006l ) 
may not apply well to high-SFR objects. Such objects, 
which presumably harbor large amounts of cold gas, may 
have higher-than-average va lues of the column density 
(e.g. JAlexander et al. 2005a). Therefore, we may under- 
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estimate the AGN fraction in such sources. 

To assess how sensitive our results are to changes in 
the population of Compton-thick sources, we altered the 
JVh distribution described in § E] by maintaining a fiat 
distribution out to TVjj = 10 25 cm~ 2 , effectively increas- 
ing the number of Compton-thick sources by a factor of 
w 3.5 (such an in crease is within the range of some recent 
predictions; e.g., iFiore et al.l 120091) . With this change, 
we find that the AGN fraction increases at all SFRs by 
~ 30%, with a somewhat larger increase at low SFRs 
or luminosities due to the larger effect of obscuration on 
the observed-frame X-ray fluxes at low redshifts (which 
generally correspond to lower SFRs and luminosities). 
Therefore, depending on the true number of Compton- 
thick sources, it is conceivable that our determination of 
the AGN fraction could be low by ~ 30% or more. 

A related issue is that of incompleteness due to statis- 
tical effects that results in underestimates of the number 
of sources near the flux limit. As noted in Section^ sim- 
ulations indicate that the incompleteness of our 70 //m 
sample is small at S/N — 3, our adopted cutoff. As a 
test of this assertion, we performed the entire analysis 
again, but with a cutoff S/N — 5. With this cutoff, 
our sample becomes almost a factor of two smaller (556 
70um sources in total versus 1022 sources for S/N = 3), 
but the trends we observe in the AGN fraction do not 
change appreciably, except for the trend with the SFR 
for higher-luminosity AGNs (io.5-8feeV > 10 43 erg s -1 ); 
see Figure [12l right). In this case, the overall trend is 
very similar, but the overall normalization is lower by 
20-30% (although they are still consistent with the 1-a 
errors). Since these higher- luminosity AGNs tend to lie 
at higher redshift, the removal of the fainter sources has 
a larger effect on this sample than on the sample of lower- 
luminosity AGNs (those with £ .5-8fceV > 10 41 erg s _1 ), 
as higher-redshift sources are fainter on average. There- 
fore, evolutionary effects could explain the lower normal- 
ization if the AGN fraction at a given SFR has a depen- 
dence on redshift. Alternatively, systematic errors in the 
estimated AGN contribution that lead to misestimates 
of the AGN contribution at shorter rest wavelengths or 
higher AGN luminosities (both of which are relevant for 



the high-redshift sources removed) could contribute to 
the observed effect. Unfortunately our data do not al- 
low us to distinguish between these two possibilities, but 
we note that longer-wavelength far-infrared observations, 
which would suffer less from AGN contamination even at 
high redshifts, would be useful in this regard. 

When we examine the growth rates of the SMBHs and 
galaxies in our sample, we find that the ratio of SMBH 
growth rate to bulge growth rate for the systems with 
identified AGNs agrees well on average with that ex- 
pected from the local scaling relation. However, when we 
include our estimates of the AGN duty cycle, the average 
SMBH growth rate is a factor of 3-10 lower (neglecting 
systems with low SFRs where incompleteness likely bi- 
ases our result significantly), suggesting that the bulges 
could be growing faster relative to the SMBHs than ex- 
pected from the local scaling. While this average SMBH 
growth rate should be considered a lower limit for the 
reasons outlined ab ove, our results are consistent with 
iMerloni et al.l (|2006f ). who derived estimates of the evo- 
lution of the total star formation and accretion density 
of the Universe. Merloni et al. found that the accretion 
rate of SMBHs falls more quickly than the SFR from 
z *~ 1 to the present day. This finding implies that the 
SMBHs were more massive relative to their bulges at this 
and higher redshifts than they are in the local universe. 
Therefore, to reproduce the local Mbuige/^BH relation, 
the bulges must grow faster relative to their SMBHs in 
the distant universe. Our estimates of the ratio of SMBH 
to bulge growth rate, corrected for the AGN duty cycle 
at a given SFR, are consistent with this scenario. 

Lastly, although a direct connection between AGN ac- 
tivity and star formation is appealing, the SFR may not 
be the fundamental quantity that drives the AGN frac- 
tion. In particular, the stellar mass of a galaxy is known 
to be related to the likelihood of AGN activity, with more 
massive galaxies on average being more likely to host an 
AGN. Additionally, simple evolution of the AGN fraction 
with redshift in all galaxies could account for our results, 
since our sample is largely degenerate between redshift 
and SFR. Recently iBrusa et all (|2009t) studied the de- 
pendence of the AGN fraction on the stellar mass for a 
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sample of X-ray selected AGNs in the CDF-S and found 
that, over a redshift range similar to that of our study, 
the AGN fraction increases from a few percent at masses 
of ~ 10 10 Mq to - 30% at masses of > 3 x 10 11 M Q . 
This trend is roughly consistent with the trend we see 
with SFR and, consequently, either property could be the 
fundamental driv er of the AGN fraction. Additionally, 
Kue et all (|2010f) found evidence that the AGN fraction 
shows a strong dependence on the stellar mass, as well as 
some dependence on redshift at a given SFR. Unfortu- 
nately, due to the difficulties in the derivation of stellar 
masses for luminous AGN hosts (due to the presence of 
significant, sometimes even dominant, emission from the 
AGN at the optical and near-infrared wavelengths typi- 
cally used for stellar-mass estimation) and our relatively 
small and flux-limited sample (which results in a degen- 
eracy of redshift and SFR), a detailed investigation of 
these questions is beyond the scope of this paper. 

In summary, given the various uncertainties discussed 
above, our findings on the total AGN fraction in distant 
luminous starbursts are broadly consistent with studies 
that have examined the ULIRG regime locally and the 
SMG regime at z ~ 2. Our results indicate, when the 
SFR is assumed to be the primary driver of AGN activity, 
that the fraction of systems hosting a higher-luminosity 
AGN increases strongly with the SFR and that SMBHs 
and bulges grow together on average, implying an inti- 
mate connection between SMBH accretion and star for- 
mation during periods of vigorous growth. 

7. CONCLUSIONS AND FUTURE WORK 

We present an analysis of AGN activity in a sample of 
starburst galaxies in the E-CDF-S and EGS fields. The 
sample was constructed from the FIDEL mid-infrared 
Spitzer survey that traces the reprocessed emission from 
young stars and AGNs. Our sample is roughly complete 
for SFRs above ~ 100 M Q yr" 1 (SFR - 600 M© yr" 1 ) to 
z = 1.0 (z = 2.0). Using this sample, we investigate how 
the incidence and strength of AGN activity, identified 
using the deep Chandra data for these fields, relates to 
the physical properties of the galaxies. In particular, 
we find that the fraction of sources that host an AGN 
depends strongly on the source's mid-infrared color, rest- 
frame mid-infrared luminosity, and, especially at higher 
intrinsic AGN luminosities, on the SFR of the galaxy. 

The dependence of the AGN fraction on the mid- 
infrared color confirms that the AGNs contribute sig- 
nificantly to the heating of the dust that emits in the 
mid-infrared. At warmer colors, approximately above 
the fiducial value o f \og(F 2 4^ m /F 70 ^ m ) > -0.7 (e.g., 
iMullanev eta l. 2009), at least - 70-80% of such sources 
host an AGN in the distant universe. The fraction of 
sources with an AGN decreases steadily with decreas- 
ing dust temperature, an indication that systems with 
cooler dust generally lack a luminous AGN and there- 
fore their mid-infrared emission is due primarily to re- 
processed emission from young stars. 

A related effect is the rise in the AGN fraction with the 
rest-frame mid-infrared luminosity, which likely suggests 
again that the AGN contributes a greater fraction of the 
emission that heats the dust at higher mid-infrared lumi- 
nosities. At the highest rest-frame mid-infrared luminosi- 
ties sampled by our study (corresponding approximately 
to ULIRG luminosities, although the selection is some- 



what different), the AGN fraction rises to ~ 60-80%, in 
rough agreement with the AGN fraction determined for 
local ULIRGs by a n umber of recent studies using a vari- 
ety of methods (e.g.. iyeilleuxll2006HImanishi et al.|[2"007t 
I Alexander et al.H2008alf ~ 

Star formation also often plays an important role in 
powering the mid-infrared emission. Therefore, we have 
estimated the contribution from the AGN to the mid- 
infrared luminosity using empirical AGN SEDs, thus al- 
lowing us to disentangle the relative contributions from 
the AGN and star formation to this luminosity. After 
accounting for the AGN-powered emission and filtering 
our sample to avoid biases, we find that star formation 
appears to power the bulk of the mid-infrared emission 
in the remaining sample and that luminous AGN activity 
is more common in systems with higher SFRs (with the 
caveat that, due to the nature of our sample, the SFR, 
mass, and redshift of our sources are largely degenerate). 
At the highest SFRs (~ 1000 M yr" 1 ), the fraction of 
sources with an AGN rises to ~ 30-40%. This fraction 
is roug hly consistent with that derived for high-re dshift 
SMGs (I Alexander et al J 120055 iLarrTetall f2010h . At 
lower SFRs (~ 30 M yr -1 ), the fraction of sources with 
a luminous AGN (£0.5-8.0 kcV > 10 43 erg s -1 ) falls to 
a few percent. However, when lower-luminosity AGNs 
(£0.5-8.0 kcV > 10 41 erg s -1 ) are included, the fraction 
is « 10% for SFR < 100 M yr -1 . This detailed re- 
lation between the AGN fraction, and hence AGN duty 
cycle, and the SFR should provide useful constraints on 
large-scale models of galaxy and SMBH ev olution (e.g., 
iDi Matteo et al.ll2008t iHopkins et al.ll2008t) . Such mod- 
els should reproduce AGN duty cycles that are consistent 
with the results presented here. 

Lastly, we made rough estimates of the growth rates of 
the SMBHs and bulges in our sample. We found that, for 
systems with detected AGN activity, the median ratio of 
bulge to SMBH growth is consistent with that expected 
from the local scaling relation, although with large scat- 
ter. This result implies that the SMBHs and bulges in 
these systems are growing concurrently on average, even 
during periods of intense star formation, at relative rates 
that would produce or maintain the scaling observed lo- 
cally. However, we do find a large scatter in this ratio, 
suggesting that in individual systems there are periods 
of rapid SMBH growth that are unaccompanied by rapid 
bulge growth (and vice- versa), although systematic un- 
certainties may account for much of this scatter. When 
the AGN duty cycle is included, the lower limit on the av- 
erage ratio of SMBH-to-bulge growth across all systems 
(not only those with detected AGN activity) suggests 
that the bulges in these distant luminous starbursts could 
be growing more quickly relative to their SMBHs than 
expected from the local scaling relation, consistent with 
recent predictions of the evolution of the to tal star forma- 
tion a nd accretion density of the Universe (|Merloni et al.l 
I2006T) . 

In summary, our results demonstrate a close connec- 
tion between AGN activity and star formation in dis- 
tant starbursts and suggest that SMBHs and their bulges 
grow together on average over a wide range of growth 
rate. However, a great deal of further work is required 
to address a number of remaining issues. For instance, 
much uncertainty remains in the determination of the 
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AGN contribution to the mid-infrared flux, a critical 
step in estimating the SFRs. Deep, longer-wavelength 
data (up to ~ 500/im) from the Herschel telescope, 
which will primarily trace cool dust emission not asso- 
ciated with the AGN, should prove invaluable in both 
increasing sample sizes and accurately estimating SFRs 
of sources at high redshift. Near- future hard X-ray ob- 
servatories, such as NuSTAR, will be helpful for the iden- 
tification of highly obscured AGNs. Lastly, improved un- 
derstanding of the relation between the observed X-ray 
emission and the bolometr ic luminosity of AGNs (e.g., 
iVasudevan fc Fabianl l200"7t ) will also be very helpful in 
constraining the AGN mid-infrared contribution. 
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APPENDIX 

PHOTOMETRIC REDSHIFTS FOR SOURCES IN THE E-CDF-S AND CDF-N 

A number of photometric redshift catalogs are available for the E-CDF-S (e.g., IWolf et al.1120041 ; IGrazian et ai]|2006D 
that have proved very useful to studies such as ours of the average properties of large samples of galaxies. Recently, 
however, new optical, UV, and near-to-mid-infrared data have become available for the E-CDF-S, as well as the 
CDF-N. These data should allow derivation of photometric redshifts that are generally of improved quality and that 
additionally include sources with fainter fluxe s. To supplement the spectroscopic redshifts used in the work described 
in this paper and others (e.g.. lXue et al.|[2010l) . we have produced photometric redshift catalogs for nearly all detected 
optical sources in the 2 Ms Chandra Deep Fields (the CDF-S and CDF-N) and the 250 ks E-CDF-S. Although we 
do not use photometric redshifts for the CDF-N in this paper, we include them here for completeness and ease of 
reference. In this appendix, we briefly describe our method of deriving photometric redshifts and present estimates of 
their quality by comparing to spectroscopic redshifts in these fields (further details are given with the catalog) Fl 

Photometric Catalogs 

For the E-C DF-S, we constructed photometric catalogs using the following catalogs: the MUSYC BVR-detected 
optical c atalog dGawiser et al.l [2006h. the COMBO-17 optical catalo g dWolf et al : l|200l l200l . the GOODS-S MUSIC 
catalog (iGrazian et all 12006ft . the MUSYC near- infrared catalog (jTavlor et all 120091 ) . the SIMPLE Spitzer IRAC 
catalog (|Damen et al.ll2009D . th e GALEX UV cata log (NUV and FUV) from the GALEX Data Release 4, and the 
GOODS-S deep J7-band catalog (jNonino et a l.l l2009l) In the CDF-N the f ollowing catalogs were used: the GOODS-N 
HST ACS and Spitzer IRAC photometric catalogs (|Dickinson et al.ll2003l) . the CDF-N Spitzer IRAC catalog derived 
from unpublished IRAC archival data, the GALEX HDF-N d eep imaging survey catalog from the GALEX Data 
Release 4, and the ACS GOODS-N region K s (< 24.5) catalog (jBarger et al.l l2008). The sources were cross matched 
by position using a matching radius of O'.'5-l" (depending on the positional uncertainty of the catalogs). The final 
E-CDF-S photometric catalog comprises a total of 105,825 unique sources, and the CDF-N catalog comprises 48,858 
sources. 

Galaxy, Hybrid, and Stellar Templates 

To model the galaxies, the 259 PEGASE galaxy templates used bv IGrazian et al. (2006)0 that span a wide range 
of star formation history and intrinsic extinction were used. To model quasars and galaxy-AGN hybrids, the 10 
AGN templates of Pollctta et al.1 ([2007)1^1 that include a variety of empirical quasar and Seyfert templates were used. 
Additionally, we constructed a set of AGN-galaxy hybrid templates by combining 5 typical AGN templates (type-1, 
type-2, QSO 1, and QSO 2) and 16 typical galaxy templates (elliptical, spiral, and starburst) with a variety of relative 
contributions. For each pair of AGN and galaxy templates, the templates were normalized by the total integrated flux 
and four hybrid templates with varying AGN/Galaxy ratios (90:10, 75:25, 50:50, and 25:75) were produced. A total 
of 330 AGN/galaxy- hybrid templates were used (330 = 5x16x4+ 10). 

In addition to the galaxy and AGN templates, 235 stellar templates from t he LePhare di stributiorPI were used to 
identify likel y stars. The templa tes used include the set of templates fr om thelPicklesI (11998T) library: the white dwarf 
templates of iBohlin et al.l {l995i); and the low- mass stellar templates of lChabrier et al.l (|2000). In addition, a set of 5 
interpolated templates was created between each of the original stellar templates, to produce a final set of 1405 stellar 
templates that should represent the full range of likely SEDs. 

Photometric Redshifts 

The Zurich Extragalactic Bayesian Redshift Analyzer (ZEBRA; Feldm ann et al.ll2006T ) was used for the photometric 
redshift derivation. Default values were used for most parameters. The reader is referred to documentation included 

17 The same method was also used to derive the photometric , . , ■ , , - c , ,. , c , , ,., 

1 1 -r. r 1, T^nc , . ,, . , . , , , , our photometric redshift estimates are ot comparable quality, 

redshifts tor the EGb used m this paper, but we have not created a 1R „ . . . .,. _ . .. . .. .. „ 

catalog of all optical sources. The reader is referred to El]for a dis- Provided with the EAZY distribution. See 

cussion of the quality of the photometric redshifts (including those , J : //m ' astr ° ■ yale • edu / eaz y/ . 

from the ECS) for the 70 fim sample presented in this paper. We See http: //cass .ucsd. edu/SWIRE/mcp/templates/swire_templates .html 

note, however, that a comparison of our EG S photometric redshifts See |http : //vwa. oamp . f r /people /arnout s/LE_PHARE . html | 
with the photometric redshifts derived by Ramos Almeida ct al. 
(2009) for 96 X-ray and mid-infrared identified AGNs shows that 
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TABLE 2 

Photometric redshift catalog for the E-CDF-S. 



RA 


Dec 


•^phot 


2 low 68% 
phot 


up 68% 
phot 


..low 95% 
phot 


up 95% 
phot 


•2spcc 


Ref. 


Template 


X-ray ID 


GOODS flag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


53.0203740 


-27.7246863 


0.105 


0.105 


0.123 


0.105 


0.172 


-1.0 


-1 


Galaxy 


-1 





53.0203700 


-27.7045750 


0.728 


0.715 


0.732 


0.694 


0.740 


0.735 


17 


Galaxy 


-1 





53.0203695 


-27.5503508 


1.271 


0.768 


1.558 


0.187 


2.525 


-1.0 


-1 


Galaxy 


-1 





53.0203590 


-27.7484474 


0.956 


0.899 


0.996 


0.858 


1.064 


-1.0 


-1 


Galaxy 


-1 






Re fere nces. — (1) IVanzella et aUflOOST: (2) iLe Fevre et all 120041; (31 fEazokolv et al.l 120041; f4) ICroom et all 1200 ll ; (5) IDickinson et all 
120041 ; (6) Ivan der Wei et aU 1 20041; (7) | Bunker et al .1120031; f8) IStanwav et afll2004li (9) IMignoli et aUI2005|; f 101 Silverm an. Mainieri, et al 
in pr epara tion; (11) ICristiani et al.||2000| ; (12) IStrolger et al.||2004t (13) IRavikumar et al.||2007l ; (14) IStanwav et al.|[2004| ; (15) ITreister et all 
1 20091; (16)|Popesso et al.||2009| (VIMOS VL T low-resolution survey); (17) |Popesso et al.||20Qg| (VIMOS VLT medium-resolution survey); (18) 
iGrazian et alj|20061 ; (191 Izheng et al.||200l 



Note. — Tablc[2]is presented in its entirety in the electronic edition. An abbreviated version of the table is shown here for guidance as to 
its form and content. The full table contains 96 columns as follows. Cols. (1-2): Source position in degrees, Col. (3): photometric redshift, 
Cols. (4-7): estimate of the 68% and 95% confidence intervals of the photometric redshift, Col. (8): spectroscopic redshift (if available), 
Col. (9): source of the spectroscopic redshift (numbers correspond to those given in the references for thi s Table), Col. (1 0): type of the 
best-fit template, Col. (11): ID of the associated X-ray source (if any) from the 2 Ms CDF-S catalog of |Luo et all 120081) or the 250 ks 
E-CDF-S catalog of ILehmer et al.l 1 120051) . Col. (12): flag indicating whether the source is inside the GOODS-S region, Cols. (13-96): the 
photometry used in the fit. 

with the catalog for futher details. In the E-CDF-S, photometric redshifts were obtained for a total of 100,318 sources 
(5507 of the original 105,825 sources had detections in fewer than 3 optical bands and were not fit). Of these, 1957 
are identified as stars, either photometrically or spcctroscopically (including 26 white dwarfs). The remaining 97,712 
sources and 649 X-ray AGNs were fit best by galaxy and AGN/galaxy hybrid templates, respectively. In the CDF-N, 
47,224 sources and 308 X-ray AGNs were fit best by the galaxy and AGN/galaxy hybrid templates, and 1323 were 
identified as stars (including 6 white dwarfs). Tables [2] and [3] give the derived photometric redshifts, the available 
spectroscopic redshifts, and the photometry used by ZEBRA for these sources. 

To assess the quality of the photometric redshifts derived by ZEBRA, comparisons are made to secure spec- 
troscopic redshifts (given Tables [2] and [3]). We used a number of quantities to assess the quality o f the photo- 
metric redshifts derived by ZEBRA: the normalized median absolute deviation (iMaronna etaDHool), CNMAD = 
1.48 x median ( | Az — median( Az)|/[l + z S peJ) , whic h gives an indication of the quality of the photometric redshifts 
after the exclusion of outliers (Brammcr et al.ll2008l ); the average absolute scatter, AAS — mean (| Az|/[1 + z spec ]) , 

which includes the effects of outliers; and the percentage of outliers with -A- — — > 0.1 and t-t — — > 0.2, where 

Az = Zphot — Zspoc- Table H] gives the quantities defined above for the photometric redshifts derived by ZEBRA for a 
number of subsamples. 

We note that, although the above indicators are commonly used to assess the quality of photometric redshifts, the 
implicit assumption in their interpretation is that the spectroscopic subsample is representative of the full sample. This 
assumption is unlikely to be entirely true, particularly when the spectroscopic sample is small relative to the number 
of sources in the total sample or, as is often the case, is brighter on average than the total sample. Additionally, 
the template-improvement step used in our derivation introduces a bias, as we have optimized the templates for the 
spectroscopic subsample. The spectroscopic subsample will therefore likely have significantly better quality than the 
full sample (unless, again, the spectroscopic subsample is fully representative). To assess the importance of these effects, 
we carried out "blind" tests for each of the three E-CDF-S subsamples (the bright galaxy sample, the faint galaxy 
sample, and the X-ray AGN sample) as follows. For each subsample, we randomly used « 3/4 of the spectroscopic 
sources for the training procedure described above and used the remaining w 1/4 of the spectroscopic sources to test 
the quality of the resulting photometric redshifts. This process was repeated eight times for each subsample to ensure 
statistically meaningful source numbers for the test sample. The results of these blind tests (see Table 0} give the 
fairest assessment of the overall quality of the photometric redshifts. In general, it appears that the use of fully trained 
subsamples gives values for AAS, cxnmad, and outlier fractions that are biased low by a factor of ~ 2-3. 

Lastly, a number of other photometric redshift catalogs exist for the E-CDF-S and CDF-N. In Table IU we compare 
the qualit y of our redshifts t o tho se from t wo widely use d phot ometric redshift catalogs: the COMBO-17 catalog of the 
E-CDF-S (|Wolf et al.l l2004. 2008) and the lCapak et~aTl ()2004h catalog of the CDF-N. It is clear that the photometric 
redshifts derived by ZEBRA are comparable or superior to those of the COMBO-17 and Capak et al. catalogs in 
most respects. We also note that the catalog of Capak et al. has an unusual deficit of sources with 2 < z < 3; 
such a deficit is not present in our ca talog. Lastly, we note that the photometric redshift catalog of E-CDF-S X-ray 
sources produced by ILuo et ail (|2010l) , which uses upper limits and deblended photometry when deriving photometric 
redshifts, supersedes the catalog presented here for E-CDF-S X-ray sources. 
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TABLE 3 

Photometric redshift catalog for the CDF-N. 



RA 
(1) 


Dec 

(2) 


-2-phot 

(3) 


low 68% 
phot 

(4) 


up 68% 
phot 

(5) 


low 95% 
phot 

(6) 


up 95% 
phot 

(7) 


-Sspcc 

(8) 


Ref. 
(9) 


Template 
(10) 


X-ray ID 
(11) 


GOODS flag 
(12) 


189.3129730 


62.3347588 


2.521 


2.406 


2.631 


2.296 


2.718 


-1.0 


-1 


Galaxy 


-1 


1 


189.3709259 


62.3344383 


0.277 


0.227 


0.326 


0.182 


0.436 


0.277 


1 


Galaxy 


-1 


1 


189.4083252 


62.3435631 


0.072 


0.020 


0.111 


0.020 


0.138 


-1.0 


-1 


Galaxy 


-1 


1 


189.3067932 


62.3343468 


0.258 


0.225 


0.304 


0.198 


0.539 


0.278 


2 


Galaxy 


-1 


1 



References — (1) IBarger et aU 120081 (2 ) [Cowie et al .1120041 : (31 IWirth et al.| [20041 ; (41 IReddy et al.|[200^ ; (51 IBarger et al.|[2003l ; (6) 
ITrouille etU]|2008t (71 IChapman et alj|2505l 



Note. — Tabic [3] is presented in its entirety in the electronic edition. An abbreviated version of the table is shown here for guidance as to 
its form and content. The full table contains 47 columns as follows. Cols. (1-2): Source position in degrees, Col. (3): photometric redshift. 
Cols. (4-7): estimate of the 68% and 95% confidence intervals of the photometric redshift, Col. (8): spectroscopic redshift (if available). 
Col. (9): source of the spectroscopic redshift (numbers correspond to those given in the references for this Table), Col. (10): type of the 
best-fit template, Col. (11): ID of the associated X-ray source (if any) from the 2 Ms CDF-N catalog of I Alexander ct al. (2003), Col. (12): 
flag indicating whether the source is inside the GOODS-N region, Cols. (13—47): the photometry used in the fit. 



TABLE 4 

Photometric redshift quality estimators. 



Case Sample # <t nm ad AAS ' Az > 0.2 ' Az| > 0.1 

iT^spcc -I 1 -''spec 



E-CDF-S field 



(a) 


all sources 


2304 


0.0116 


0.0262 


2.00% 


3.04% 


(b) 


thr < 24 galaxies (trained) 


1699 


0.0118 


0.0215 


1.12% 


1.88% 


(c) 


m R < 24 galaxies (blind) 


1738 


0.0345 


0.0501 


3.86% 


8.46% 


(d) 


mR < 24 galaxies (blind-C17) a 


1530 


0.0291 [0.0244] 


0.0404 [0.0574] 


2.87% [7.12%] 


6.00% [9.61%] 


(c) 


thr > 24 galaxies (trained) 


605 


0.0113 


0.0396 


4.46% 


6.28% 


(f) 


mR > 24 galaxies (blind) 


619 


0.0612 


0.1142 


14.05% 


25.20% 


IK) 


m R > 24 galaxies (blind-C17) a 


230 


0.0500 [0.1057] 


0.0607 [0.1811] 


3.81% [25.22%] 


14.62% [40.43%] 


(h) 
(i) 


X-ray AGNs (trained) 


283 


0.0094 


0.0140 


1.06% 


2.12% 


X-ray AGNs (blind) 


315 


0.0436 


0.0931 


14.92% 


24.44% 


(j) 


X-ray AGNs (blind-C17) a 


217 


0.0495 [0.0251] 


0.0835 [0.0738] 


13.52% [11.98%] 


22.53% [20.28%] 


CDF-N field (GOODS-N region) b 


(k) 


all sources 


2672 


0.0229 [0.0440] 


0.0480 [0.0819] 


4.27% [8.42%] 


8.68% [16.84%] 


(1) 


tur < 24 galaxies 


1837 


0.0227 [0.0362] 


0.0409 [0.0497] 


2.78% [3.54%] 


6.26% [9.80%] 


(m) 


run > 24 galaxies 


835 


0.0234 [0.0753] 


0.0637 [0.1526] 


7.54% [19.16%] 


14.01% [32.34%] 


(n) 


X-ray AGNs 


164 


0.0142 [0.0760] 


0.0380 [0.1505] 


3.66% [20.73%] 


7.32% [26.22%] 



CDF-N field (non-GOODS-N region) b 



(o) 


all sources 


2687 0.0245 


0.0440 


0.0538 


0.0819 


5.02% 


8.41%] 


10.98% [16.86%] 


(p) 


mR < 24 galaxies 


1848 0.0256 


0.0364 


0.0478 


0.0497 


3.63% 


3.52%] 


9.58% [9.85%] 


1' 


ran > 24 galaxies 


839 0.0226 


0.0744 


0.0671 


0.1528 


8.10% 


19.19%] 


14.06% [32.30%] 


(r) 


X-ray AGNs 


212 0.0118 


0.0804 


0.0457 


0.1617 


4.72% 


21.70%] 


7.08% [29.72%] 



a To obtain these numbers, we used sources with spectroscopic redshifts that have photometric rcdshifts in both the ZEBRA and 
COMBO-17 catalogs. Values in brackets arc derived from the COMBO-17 catalog. 

k Values given for the CDF-N arc obtained from the fully trained subsamplcs. Values in brackets are derived from the photometric 
redshift catalog of Capak ct al. (private communication) based on the photometric catalog of Capak ct aQ J 200411 - 



